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Preamble 

As in earlier reports, we will continue to break our effort into six distinct units: 

• Atmospheric Correction Algorithm Development 

• Whitecap Correction Algorithm 

• In-water Radiance Distribution 

• Residual Instrument Polarization 

• Pre-launch/Post-launch Atmospheric Correction Validation 

• Detached Coccolith Algorithm and Post-launch Studies 

This separation has been logical thus far; however, as launch of AM-1 approaches, it must be recognized 
that many of these activities will shift emphasis from algorithm development to validation. For example, the 
second, third, and fifth bullets will become almost totally validation-focussed activities in the post-launch 
era, providing the core of our experimental validation effort. Work under the first bullet will continue into the 
post-launch time frame, driven in part by algorithm deficiencies revealed as a result of validation activities. 

Prior to the start of the 1999 fiscal year (FY99) we were requested to prepare a brief plan for our FY99 
activities. This plan is included as Appendix I. The present report describes the progress made on our 
planned activities. 
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Abstract 

Significant accomplishments made during the present reporting period: 

• Submitted version 4 of the Water-leaving Radiance ATBD. 

• Submitted version 4 of the Detached Coccolith ATBD. 

• Completed analysis of a whitecap data set acquired in the Tropical Pacific. The results 
showed (1) the whitecap augmented reflectance in the visible for a given wind speed in 
the 8-12 m/s range was about 2.5 times less than that predicted by whitecap reflectance 
models, and (2) the augmented reflectance in the NIR (860 nm) was approximately 80% 
that in the visible, and only slightly dependent on the wind speed. 

• Participated in pre-INDOEX and INDOEX field experiments. Obtained sun photometer 
and micro pulse lidar data along 55 days of ship tracks, and in- water up welling radiance 
distribution data at 27 stations. 
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Atmospheric Correction Algorithm Development 


As we have successfully processed SeaWiFS imagery (reformatted to MODIS) using our MODIS code, we 
can be assured that the implementation is correct. Thus, we focussed on ensuring the correctness of the basic 
correction algorithm, improving the basic algorithm to provide for an accurate correction in the presence of 
strongly-absorbing aerosols, and addressing several MODIS-specific issues. 


1. Validate /Imvrove the Basic MODIS Algorithm using SeaWiFS 

We will begin the algorithm validation and improvement process using SeaWiFS imagery. 

a. Task Progress: 

We have been processing SeaWiFS imagery using both the SeaDAS processing software and the University 
of Miami processing system DSP (which affords us more control over the processing). The main focus of 
this processing has been the development of newer versions of the algorithm, and this is discussed in a later 
section. 

We have processed significant amounts of imagery that show that the “6-8” band selection for SeaWiFS 
atmospheric correction retrieves water-leaving radiances that are less noisy than the “7-8” band selection. 
There are several reasons why this should be true (e.g., the 7-8 algorithm requires a correction to the band 
7 radiances for O 2 absorption near 759 nm, and bands 7 and 8 are closer together than bands 6 and 8 which 
decreases the influence of sensor noise on the choice of candidate aerosol models); however, we believe that 
use of the 7-8 algorithm should be continued because there is less contaminating water-leaving radiance in 
band 7 than in band 6. This is relevant to the MODIS algorithm in that MODIS uses bands 15 and 16, 
which are similar to the SeaWiFS bands 7 and 8 except that there is no O 2 correction required for MODIS. 

Recently, in collaboration with R. Evans, we have started to participate in a SeaWiFS effort to understand 
the causes of the negative water-leaving radiances in the blue (412, 443, and 490 nm) retrieved under certain 
conditions. Although relatively few images have been examined, it appears that this problem is associated 
with one of several occurrences: (1) a very clear atmosphere; (2) absorbing aerosols; (3) Case 2 waters; and 
(4) very large scan (viewing) angles. 

In the case of item (1) we have found instances in which the residual radiance after removal of the Rayleigh 
component was very small, or even negative, (i.e., p*(865) — p r (865) < 0). The most likely cause of this 
behavior is a radiometric calibration in band 8 (865 nm). All of the radiometric calibration efforts relating 
to SeaWiFS to date have assumed that there is no error in the calibration of band 8, e.g., Gordon et al 
[1998]. This assumption means that bands 1 through 7 are calibrated relative to band 8. However, it is 
clear that a calibration error of a few % could exist in band 8. Thus, we arbitrarily assumed that SeaWiFS 
band 8 was calibrated 5% too low, and used the surface and atmosphere data acquired during the SeaWiFS 
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initialization cruise to “recalibrate” the sensor. This effort resulted in fewer negative water-leaving radiances 
and eliminated most instances in which p*(86 5) - p r (86 5) < 0. 

In the case of item (2), we are developing new versions of the correction algorithm that will address absorbing 
aerosols. These are discussed in a later section. 

In the case of items (3) and (4), the atmospheric correction algorithm is being applied in situations in which 
it was not designed to operated. For Case 2 waters, the “spectral optimization algorithm” described in a 
later section my provide a processing alternative. The large scan angle problem (item (4)) is not an issue 
with MODIS because of a more restricted scan and the absence of sensor tilt. 

b. Anticipated Future Actions: 

We will continue to test the algorithm with SeaWiFS imagery. In particular, we will continue to examine 
the problem of negative water-leaving radiances, and will continue work on the “spectral matching” and 
“spectral optimization” algorithms as described below. 

c. Publications: None. 

2. Implement and Test the “Spectral Matching Algorithm ” 

We will continue our basic study of the “spectral matching ” algorithm. Unlike the basic MODIS algorithm , the 
spectral matching algorithm has the potential for atmospheric correction in the presence of strongly- absorbing 
aerosols. 

a. Task Progress: 

Gordon [1997] has described the difficulties of atmospheric correction in the presence of absorbing aerosols. 
These include the fact that (1) it is impossible to distinguish between weakly and strongly absorbing aerosol 
using NIR bands because the spectral variation of the aerosol reflectance there depends mostly on the 
aerosol’s size distribution, (2) some aerosols, e.g., Saharan dust are nearly nonabsorbing in the green, red, 
and NIR, but absorb strongly in the blue, (3) the impact of aerosol absorption depends strongly on the 
aerosol’s vertical distribution, and (4) the impact of aerosol absorption is most significant in the blue and 
blue-green bands where Rayleigh scattering is largest. These facts imply that aerosol absorption can only be 
ascertained through measurements in the blue region of the spectrum — precisely the region most sensitive 
to the water’s chlorophyll concentration. Thus, atmospheric correction cannot be effected in the absence 
of knowledge regarding the water-leaving radiance throughout the spectrum, which requires a model of the 
marine reflectance as a function of the chlorophyll concentration. Such models exist (e.g., Gordon et al. 
[1988]), and have been incorporated into two atmospheric correction procedures that we refer to as the 
“spectral matching algorithm” and the “spectral optimization algorithm.” 

The spectral matching algorithm (SMA) [Gordon, Du and Zhang , 1997] and the spectral optimization al- 
gorithm (SOM) [Chomko and Gordon , 1998] have been implemented within the SeaDAS image processing 
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environment. Both algorithms perform atmospheric correction and retrieve water properties simultaneously, 
and represent a significant departure from the basic MODIS algorithm. The SOA has been tested more 
throughly, because we believe that it is a more versatile in the sense that it does not require accurate aerosol 
models to effect a correction. However, in the long run the SMA will likely be more useful in that, given 
realistic aerosol models, it will allow an accurate estimate of the aerosol optical thickness. 

We tested the SOA using SeaWiFS imagery and compared its performance with the basic MODIS algorithm. 
The results are presented in detail in Appendix II. Briefly, the imagery used the performance of the 
algorithms on a turbid day followed (two days later) by a very clear day in the Middle Atlantic Bight 
(MAB). As the water properties were not expected to change significantly over the two-day period of the 
imagery (October 1998), the stability of the derived pigment concentrations over the period provided an 
excellent test of the algorithms. It was found that the SOA and the MODIS algorithms performed equally 
well on the clear day, but the SOA was considerably better on the turbid day. This leads us to believe that 
the SOA approach is a viable alternative to the standard MODIS algorithm. It also provides what may be 
a viable approach to atmospheric correction on Case 2 waters. 

b. Anticipated Future Actions: 

We will continue to test the SOA under a variety of conditions; however, we will also continue to develop 
the SMA. There are several reasons why the SMA is more desirable than the SOA: (1) the SMA is capable 
of retrieving the aerosol optical depth if the candidate aerosol models are realistic; (2) the SMA is capable 
of being incorporated into the standard MODIS algorithm (at reduced resolution) in a relatively simple 
manner; (3) there are situations (e.g., Saharan dust) where the SOA will not perform well because of the 
assumption that the index of refraction of the aerosol is independent of wavelength; (4) as realistic models 
of aerosols are developed, they can be incorporated into the SMA environment in a simple manner, and (5) 
at a minimum it will be easy to use the SMA to provide a flag within the MODIS algorithm that will signal 
the probable presence of absorbing aerosols, and indicate that the quality of the derived products cannot 
be assured. We expect the SMA to be the approach of choice for atmospheric correction in the presence of 
wind-blown dust. 

c. Publications: None. 

3. Optical Properties of Wind-Blown Dust 

We need quantitative estimates of the optical properties of wind-blown dust, e.g., from Africa, to provide a 
proper model to effect atmospheric correction using the standard MODIS algorithm. 

a. Task Progress: 

The optical properties of wind-blown dust are required to effect atmospheric correction using the standard 
MODIS algorithm, and possibly the spectral matching algorithm as well. In collaboration with post doctoral 
associate Cyril Moulin, the approach we have taken for development/improvement of a wind-blown dust 
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model is to utilize SeaWiFS imagery acquired off the African Coast (east of the Canary Islands and west 
of Dakar). In this region, the dust is mixed high in the atmosphere, and the optical thickness is often very 
large (r a (865) > 2) so the marine reflectance is small compared to the reflectance of the aerosol. In this 
situation, small errors in the marine reflectance will have an insignificant impact on the retrieved aerosol 
reflectance. We used marine reflectances computed from the Gordon et al. [1988] semi-analytic model for the 
marine reflectances corresponding to the SeaWiFS monthly mean chlorophyll concentrations, and extracted 
aerosol reflectances 

Pt( A) - Pr( A) - t(X)pt(X) = pa{ a), 

where the subscripts £, r, and w , refer to total, Rayleigh, and water, respectively, and t is the diffuse 
transmittance of the atmosphere, along a track that passed through a dust plumb of varying opacity. The 
desert dust model size distribution proposed by Shettle [1984] was combined with the imaginary part of the 
refractive index proposed by Patterson [1981] to compute the variation of pa( A) as a function of p^(865), 
for comparison with the observed values along the track. We call this the BDS model (Background Desert 
Shettle). This model was found by Moulin et al [1997] to fit sun photometer aerosol optical depth data 
better than several other popular models in the literature. We found that the fit to the spectral data was 
very poor. This led us to believe that the Patterson [1981] absorption indices are too large. We are now 
trying a variety of absorption indices and modifications to the Shettle [1984] size distribution to develop a 
better dust model. 

We tested the efficacy of a variant of the SMA for atmospheric correction in a region off the coast of Africa 
that was subjected to a high concentration of Saharan dust. In this case we used the BDS model with 
varying concentrations of large particles, and with varying thicknesses for the aerosol layer. The results were 
encouraging and suggested that atmospheric correction can be effected under dust conditions with aerosol 
optical thicknesses at 865 nm as high as unity. 

b. Anticipated Future Actions: 

We will continue to refine the dust models described above with the goal that they will be available for 
inclusion in the MODIS correction algorithm, or its SMA variant, by the end of the next reporting period. 
As described above, preliminary attempts to use the BDS model in atmospherically correcting imagery off 
the African coast have shown considerable promise. These will continue. 

We also have a unique opportunity to use SeaWiFS imagery that was acquired over the Caribbean in 1998 
during particularly intense dust periods of July and August, simultaneous with our surface measurements of 
the aerosol vertical profile (using our MPL) and aerosol optical thickness (using a CIMEL sun/sky radiome- 
ter) made from St. Johns, VI, to validate our model for African dust. 

c. Publications: None. 

1. MODIS-Svecific Issues 
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We need to address the detection and removal of thin cirrus clouds, methods for efficiently including earth- 
curvature effects, out-of-band corrections, BRDF effects on the diffuse transmittance, correct for polarization 
sensitivity of MOD IS, and the efficacy of atmospheric correction for removal of the aerosol effect from the 
measurement of the fluorescence line height. 

a. Task Progress: 

We have used the MODIS relative spectral response (RSR) functions to complete the out-of-band corrections 
to the algorithm using the Gordon [1995] methodology. These will be included in the at-launch version of 
the algorithm. 

We have included the MODIS polarization-sensitivity correction in the MODIS processing code. This task 
will be complete when the SBRS/MCST polarization-sensitivity data are received from MSCT. 

We are working to develop a model of the subsurface BRDF. These efforts are described below (under 
“In- water Radiance Distribution”). 

b. Anticipated Future Actions: 

We know that there are difficulties atmospherically correcting SeaWiFS imagery when the scan angle exceeds 
~ 50°. Although scan angles greater than 50° comprise only a minor portion of the MODIS scan, the larger 
the scan angle that can be successfully corrected, the larger the base of MODIS data that will be available. 
The primary of cause of algorithm failure at large scan angles may be the fact that earth-curvature is 
ignored in the algorithm. Thus, we will investigate the requirements of including earth-curvature effects in 
the MODIS algorithm following the Ding and Gordon [1994] approach. This involves a modification of the 
look-up tables for the top-of-the-atmosphere contribution from Rayleigh scattering. However, before actually 
embarking on an implementation, we will examine SeaWiFS imagery at high latitude and large scan angle 
to assess the impact of neglecting earth-curvature in the algorithm. 

Finally, because of uncertainty in the performance of MODIS Band 26, and because we need to assess 
whether our radiative transfer codes are sufficiently accurate to study removal of the aerosol effect from the 
measurement of the fluorescence line height, we will examine this issue and thin cirrus clouds only if time 
permits. 

c. Publications: None. 
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Whitecap Correction Algorithm 


As the basic objectives of the experimental portion of this task has been realized (acquiring whitecap radio- 
metric data at sea), experimental work is being suspended until the validation phase, except insofar as the 
radiometer is being operated at sea when it is sufficiently important to do so, e.g., it was operated during 
the SeaWiFS Initialization Cruise (MOCE- 4). Our basic goal for the rest of the project is to maintain 
experience in operating and maintaining the instrumentation in preparation for the validation phase of the 
contract, and to complete the analysis of the whitecap data acquired thus far. 

5. Maintaining Measurement Cayability 

We need to maintain our ability to make whitecap measurements. 

The instrumentation and data acquisition system has been maintained and will be ready for the MODIS 
Initialization Cruise. 

6. Reduction and Analysis of Existing Data 

We need to complete the analysis of the data already acquired. 

a. Task Progress: 

We have concentrated on a reanalysis of the Tropical Pacific whitecap data set that was acquired in 1996. 
This is a unique data set, as it was acquired in the trade winds with moderately high winds (8-12 m/s) and 
practically unlimited fetch and duration. We developed an alternative method of analysis that we believe is 
more robust and provides greater confidence in the results. The complete analysis is provided as Appendix 
III. Briefly, there were two significant results that emerged: 

1. The whitecap augmented reflectance in the visible for a given wind speed in the 8-12 m/s 
range was found to be about 2.5 times less than that predicted by whitecap reflectance 
models [Gordon and Wang , 1994; Koepke , 1984]. 

2. The augmented reflectance in the NIR (860 nm) was approximately 80% that in the 
visible, and dependent on the reflectance (i.e., whitecaps that were brighter in the visible 
showed a greater reflectance decrease in the NIR than whitecaps with low reflectance in 
the visible) . This NIR decrease is in agreement with reflectance observations of breaking 
waves in the surf zone [Frouin, Schwindling and Deschamps , 1996] and of ship-wake foam 
[Moore, Voss and Gordon , 1998]. 

Last year, the SeaWiFS Project informed us that the present whitecap algorithm (See MODIS Normalized 
Water-leaving Radiance Algorithm Theoretical Basis Document, Version 4, April 30, 1999) was causing the 
atmospheric correction algorithm to fail in the South Atlantic. We provided an algorithm adjustment (based 
on our earlier analysis of the Tropical Pacific whitecap data), and were informed by them that the algorithm 
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appeared to be working much better after the adjustment. That adjustment was similar, but not identical, 
to the one resulting from our reanalysis. 

b. Anticipated Future Actions: 

We will operate the whitecap radiometer during the MODIS Initialization and Validation Cruises. 

c. Publications: 

Spectral Reflectance of Whitecaps: Their Contribution to Water-Leaving Radiance, K.D. Moore, K.J. Voss, 
and H.R. Gordon, Jour. Geophys. Res. (Submitted). 
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In-water Radiance Distribution (BRDF) 


During FY ’99 our objectives for this task are: maintaining experience in operating and maintaining the 
instrumentation in preparation for the validation phase of the contract; acquiring more field data; and using 
the field data to develop a model of the BRDF as a function of the solar zenith angle and the water’s 
chlorophyll concentration. 

7. Maintaining Measurement Capability 

We need to maintain our ability to make BRDF measurements. 

The instrumentation and data acquisition system has been maintained (actually they have been deployed in 
two seperate field experiments) and will be ready for the MODIS Initialization Cruise. 

8. Acquisition of More Field Data 

To build and thoroughly test a BRDF model, we need to acquire data over a wider range of chlorophyll 
concentrations. 

We deployed the radiance distribution camera (RADS) during INDOEX in the Indian Ocean in February- 
April 1999. The data are presently being reduced to provide the BRDF as a function of chlorophyll concen- 
tration during the cruise. 

9. Build a BRDF Model 

We need to develop a model relating the subsurface radiance distribution (BRDF) to the chlorophyll concen- 
tration of the water. 

a. Task Progress 

Our efforts to model the BRDF have yet to yield an acceptable model. We started by using an ocean model 
developed by Morel and coworkers [Morel and Gentili , 1991; Morel and Gentili , 1993; Morel and Gentili , 
1996; Morel, Voss and Gentili , 1995]. This model, although often utilized in the literature, was not capable 
of reproducing our measurements of the distribution of subsurface upwelling radiance L u (0,cj)), normalized 
to nadir, for a chlorophyll concentration near 0.1 mg/m 3 with an accuracy sufficient to be useful in remote 
sensing. We believed that the departure of the computed L u (6,(j)) from the experimental distribution was 
due to the neglect of Raman scattering. Thus we undertook a reexamination of the contribution of Raman 
scattering to L u {6,(j)) in the light of recent measurement of the absorption coefficient of pure water [ Pope 
and Fry , 1997] and the spectral variation of the Raman-scattering cross section [Bartlett et al., 1998]. The 
results of this study [ Gordon , 1999] showed that Raman scattering made a larger contribution than expected 
on the basis of earlier computations [ Waters , 1995], and that the fraction of L u {6 , (j>) contributed by Raman 
scattering in the blue did not significantly decrease with increasing chlorophyll concentration. Unfortunately, 
accounting for the effects of Raman scattering did not significantly improve the fit of the computed L u {6 , </>) 
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to the experimental data. Therefore, the deficiency of the model must be elsewhere. There are two leading 
candidates: (1) the scattering phase function used for the particles is unrealistic; and (2) the assumption 
that scalar radiative transfer theory (polarization is ignored) is adequate is not correct. A third possibility 
is that the measured radiances are in error due to instrument self shading. 

b. Anticipated Future Actions: 

The phase function used to characterize the particle scattering was actually measured in-situ and has long 
been used to characterize particle scattering in the ocean [ Petzold , 1972]. Therefore, we are looking first 
toward the use of scalar radiative transfer theory in the computations. Although very few computations 
have been carried out in the vector mode, the errors in the radiance computed using scalar theory can 
be roughly the same magnitude as encountered in our modeling of the BRDF [Kattawar and Adams , 1989; 
Kattawar and Adams , 1990; Kattawar, Adams and Tanis , 1988]. Thus, before experimenting with the particle 
phase function, we are including polarization in our computation of both the elastic and inelastic (Raman) 
scattering components. When this is complete, we will test the model again against the data acquired during 
the SeaWiFS initialization cruise and, if successful, against the data acquired during the INDOEX cruise. 

In addition, we are collaborating with J.P. Doyle at JRC, Ispra, Italy, in modeling the self-shading charac- 
teristics of our BRDF instrument (RADS). This study uses a 3-d Monte Carlo radiative transfer code to 
compare the radiance distribution in the presence and the absence of RADS as a function of the optical 
properties of the water. 
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Pre-launch/Post-launch Atmospheric Correction Validation 


The original objectives of this task were fourfold: (1) study aerosol optical properties over the oceans to 
assess the applicability of the aerosol models proposed for use in the atmospheric correction algorithm; 
(2) measure the aerosol optical properties (including their vertical distribution) from a ship during the 
initialization/ calibration/^ validation cruises; (3) determine how accurately the radiance at the top of the 
atmosphere can be estimated based on measurements of sky radiance and aerosol optical thickness at the sea 
surface (i.e., vicarious calibration); and (4) utilize data from other sensors that have achieved orbit (OCTS, 
POLDER, SeaWiFS ...) to validate and fine-tune the correction algorithm. We have obtained a significant 
amount of data toward (1), designed, constructed, or purchased, instrumentation to acquire data for (2), 
completed (3), and (4) was discussed under the first activity above (Atmospheric Correction Algorithm 
Development). Objectives for the post-launch validation phase are: maintain experience in operating and 
maintaining the instrumentation in preparation for the validation phase of the contract; complete analysis 
of data already acquired; and participate in the validation phase. 

10. Maintaining Measurement Capability 

We need to maintain our ability to make atmospheric measurements at sea. 

a. Task Progress: 

We exercised our instrumentation - sun photometers, sky radiance camera, aureole camera, and micro pulse 
lidar (MPL) - during the INDEOX and the pre-INDOEX Cruise aboard the RV Ron Brown. Our MPL 
worked well during the pre-INDOEX cruise, but failed on the last day. As there were two other MPL units 
on the INDOEX leg, and E.J. Welton of our group was able to keep them operating, no data were lost 
because of our MPL failure. The aureole camera failed because of a disfunctioning CCD array, and was lost 
during the entire period. The sky radiance camera did not function during pre-INDOEX, was repaired for 
INDOEX, but then failed again. However, little data were lost as there was only one opportunity to acquire 
high-quality data, which requires a clear sky at large solar zenith angles (> 60°). 

Our CIMEL sun/sky scanning radiometer has been repaired, recalibrated, and reinstalled at Fort Jefferson 
in the Dry Tortugas. It is now transmitting data to the AERONET network. 

b. Anticipated Future Actions: 

The MPL has been sent back to the manufacturer for an upgrade of the optical design that will prevent 
detector damage when the laser becomes misaligned (the cause of the failure during INDOEX). The aureole 
camera is being repaired. We expect our measurement capability to be intact for the MODIS Initialization 
Cruise. 

11. Complete Analysis of Existing Data 

We need to complete the analysis of the data acquired previously. 
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a. Task Progress: 

During much of this reporting period our personnel were at sea participating in pre-INDOEX and in IN- 
DOEX. The highest priority was collecting the data during the field experiments, and the next priority was 
reduction and analysis of the data that was collected. Thus, analysis of existing data was placed on the 
“back burner” in favor of the new data sets. Because of this, the MPL data from the INDOEX cruises has 
been completely reduced along with the sun photometer data. 

b. Anticipated Future Actions: 

We will reduce the the remaining data from INDOEX next, after which we will work on data acquired in 
Hawaii during the MOCE-4 cruise, in particular the sky camera and aureole data. 

Two new graduate students have joined our group to take over the responsibilities of E.J. Welton and J. 
Ritter who have left to work at GSFC and MSFC, respectively. Graduate student David Bates has learned 
how to operate the CIMEL system and reinstalled this system in the Dry Tortugas. He is also learning the 
MPL system and will work with it when it returns from the manufacturer. Finally, he is also beginning 
analysis of the existing CIMEL data. Graduate student Hong Du is learning to operate the aureole system 
and helping in the redesign of the instrument. We are currently rebuilding this instrument to improve its 
operation in the field. Along with learning to operate the instrument he will be investigating what options 
exist for reducing this data to aerosol properties. 

c. Publications: 

E.J. Welton, K.J. Voss, D.L. Savoie, and J.M. Prospero, “Measurements of Aerosol Optical Depth over the 
North Atlantic Ocean: Correlations with Surface Aerosol Concentrations” (Under revision to Jour. Geophys. 
Res.). 

J.M. Ritter and K.J. Voss, “ A new instrument to measure the solar aureole from an unstable platform.” 
(Revised to Jour. Atmos. Ocean. Tech.). 

B. Schmid, P.B. Russell, J.M. Livingston, S. Gasso, D.A. Hegg, D.R. Collins, R.C. Flagan, J.H. Seinfeld, 
E. Ostrom, K.J. Noone, P.A. Durkee, H.H. Jonsson, E.J. Welton, K.J. Voss, H.R. Gordon, P. Formenti, 
M.O. Andreae, V.N. Kapustin, T.S. Bates, and P.K. Quinn, “Clear column closure studies of urban-marine 
and mineral-dust aerosols using aircraft, ship, and ground-based measurements in ACE-2,” (ALPS99, 18-22 
January 1999, Meribel, France). 

B. Schmid, J.M. Livingston, P.B. Russell, P.A. Durkee, H.H. Jonsson, D.R. Collins, R.C. Flagan, J.H. 
Seinfeld, S.A. Gasso, D.A. Hegg, E. Ostrom, K.J. Noone, E.J. Welton, K.J. Voss, H.R, Gordon, P. Formenti, 
and M.O. Andreae, “Clear sky closure studies of lower tropospheric aerosol and water vapor during ACE-2 
using airborne sun photometer, airborne in-situ, space-borne, and ground-based measurements.” (Submitted 
to Tellus Special Issue on ACE-2) 
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E.J. Welton, K.J. Voss, H.R. Gordon, H. Maring, A. Smirnov, B. Holben, B. Schmid, J.M. Livingston, 
P.B. Russell, P.A. Durkee, P. Formenti, and M.O. Andreae, “Ground-based Lidar Measurements of Aerosols 
During ACE-2: Instrument Description, Results, and Comparisons with other Ground-based and Airborne 
Measurements.” (Submitted to Tellus Special Issue on ACE-2) 

12. Post-launch Validation 


We will participate in MODIS post-launch validation. 

a. Task Progress: 

We used the MPL and in-water radaince distribution system (RADS) during pre-INDOEX and INDOEX 
in January- April of 1999. In this case, SeaWiFS imagery will be used as a surrogate for MODIS in the 
validation program. The INDOEX cruise (and the associated pre-INDOEX, or ACE cruise) offered a unique 
opportunity to collect aerosol optical information, satellite data (SeaWiFS), and aerosol chemistry data. The 
pre-INDOEX cruise went from Norfolk, VA to Cape Town, South Africa. During this cruise, we were able to 
sample an extensive bio-mass burning plume that extended off of Africa into the Atlantic. We obtained MPL 
data along the ship track, which allows us to look at the vertical distribution of the aerosols, an important 
piece of information for these absorbing aerosols. In addition other workers on the ship were making chemical 
and other optical measurements along this track. 

During the INDOEX cruise the focus was on the extensive pollution plume, produced on the Indian sub- 
continent, which extends over the Indian Ocean. The cruise track went from the clean air in the southern 
portion of the Indian ocean to the polluted air near India. We were able to collect the MPL information 
for this whole track along with RADS data. Other researchers on the ship measured additional in-water 
optics and aerosol chemistry. We were also obtaining SeaWiFS LAC data. For the clear days this data will 
be useful to investigate the performance of the atmospheric correction algorithms for this heavily polluted 
environment. It will also be valuable for developing models for the aerosol pollution in a manner similar to 
that described for Saharan dust in Section 3. 

We ended up with approximately 55 days, or about 16,000 nautical miles of MPL data. In addition we had 
27 in- water radiance distribution casts, in varying water types. 

b. Anticipated Future Actions: 

Given our experience with the aureole system during INDOEX we have repaired the camera portion and 
are rebuilding much of this system to improve operation. These repairs will be done in time for an expected 
MODIS initialization cruise. The other major effort will be in reducing the backlog of data we have acquired 
during the recent cruises. 

Finally we anticipate participating in the MODIS initialization cruise. The timing of this cruise will, of 
course, depend on the launch of TERRA. 
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Detached Coccolith Algorithm 


The objectives of our FY 99 effort are to complete analysis of cruise work done to date, continue new field 
data in the Gulf of Maine and Gulf of Mexico, use these data to improve algorithm performance and to 
apply and validate the coccolithophore algorithm using SeaWiFS data. 

IS. Processing of completed yre-launch cruises 

a. Task Progress: 

As of June 1999, all of our 1998 pre-launch cruise data has been processed, at least to first-order. Tables 
of time, inherent optical properties, and apparent optical properties have been submitted to the NASA 
SeaBass data base. Discrete samples from the 1998 field season have been worked up, with the exception 
of the suspended calcite data (due to lack of available time on the Graphite Furnace Atomic Absorption 
Spectrometer at the University of Maine). For this current “coccolithophore season,” we have completed 7 
transects across the Gulf of Maine. Discrete samples for this work have yet to be processed, except for the 
chlorophyll a concentrations. Another 7-day pre-launch cruise has just been completed aboard the research 
vessel Linke, in the Gulf of Maine in which coccolith concentration, water-leaving radiance, and inherent 
optical properties (including acid-labile backscattering) were measured. 

b. Anticipated Future Actions: 

Running of the suspended calcite samples from the 1998 field season will probably be done after this busy 
summer field season. Discrete samples from the 1999 field season will be processed as they become available. 

c. Publications: 

Balch, W. M., D. Drapeau, and J. Fritz, Monsoonal forcing of calcification in the Arabian Sea, Deep Sea 
Res. II. (Accepted, in revision). 

11. New field efforts 

a. Task Progress: 

We have 13 more days of shiptime across the Gulf of Maine in 1999. These will specifically involve sampling 
for coccoliths and their related backscattering. We have added a Hobi Labs Hydroscat 2 to our underway 
sampling system. After an initial failure of the instrument’s mother board, it is now functioning. This 
instrument will provide an estimate of the wavelength dependence of backscattering (at 470 and 680 nm). 

b. Anticipated Future Actions: 

We are involved in the planning of a coccolithophore cruise south of Iceland, for either summer 2000, or 
2001. Remaining discrete samples for PIC, chlorophyll, and POC will be processed. 
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c. Publications: None. 

15. Improving algorithm performance 

a. Task Progress: 

A second manuscript on our Arabian Sea results is in preparation at the time of this writing. This manuscript 
will specifically deal with the calibration of the acid-labile backscattering vs suspended calcite relationship in 
this region of known coccolithophore production and sedimentation. Moreover, the first continuous record of 
particulate inorganic carbon and particulate organic carbon, measured over ~ 3500 kilometers of ship track, 
will be presented. 

b. Anticipated Future Actions: 

Completion of manuscript, and beginning of the Gulf of Maine algorithm work. 

c. Publications: 

Balch, William M., David T. Drapeau, Terry L. Cucci, and Robert D. Vaillancourt, Katherine A. Kilpatrick, 
Jennifer J. Fritz, Optical backscattering by calcifying algae-Separating the contribution by particulate inor- 
ganic and organic carbon fractions J. Geophys. Res., 104C, 1541-1558 (1999). 

Balch, W. M., D. Drapeau, B. Bowler, and J. Fritz. Continuous measurements of calcite-dependent light 
scattering in the Arabian Sea. Deep Sea Res. II. (In Preparation). 
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Appendices 


I. OCEAN OBSERVATIONS WITH EOS/MODIS Algorithm Development and Post Launch Studies: Plans 
for FY 99. 

II. Atmospheric Correction of Ocean Color Imagery: Use of the Junge Power-Law Size Distribution. 

III. Spectral Reflectance of Whitecaps: Their Contribution to Water-Leaving Radiance. 
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Preamble 

The coming Fiscal Year (1999) was to be heavily focused on validation of MODIS-derived 
products. However, the delay of the launch of EOS AMI requires some modification of 
the plan. Our approach for the coming year is to use SeaWiFS for validating algorithms 
(rather than MODIS-derived products) in a manner similar to the way the MODIS Land and 
Atmospheres Groups have already been using the MAS. In addition, the delay in the 
launch will allow us additional time to enhance the already existing algorithms prior to 
launch. 

We break our effort for FY 99 into five distinct (although interrelated) units: 

• Atmospheric Correction Algorithm Development 

• Whitecap Correction Algorithm 

• In-water Radiance Distribution (BRDF) 

• Pre-launch/Post-launch Atmospheric Correction Validation 

• Detached Coccolith Algorithm 

In what follows, it will be seen that addressing the various tasks within these units requires 
fundamental studies, a through examination of SeaWiFS imagery, use of SeaWiFS 
imagery as a tool, and a maintenance of measurement and data-analysis capability through 
the prelaunch era into the formal validation phase. Developing the required capability of, 
and experience in, processing SeaWiFS imagery will also facilitate the efficient handling of 
MODIS imagery (e.g., for QA) in the post-launch era. 


Atmospheric Correction Algorithm Development 

Objectives and Proposed Activities (FY '99) 

During FY '99 there are several objectives under this task. They are focused on ensuring 
the correctness of the basic correction algorithm and its implementation, improving the 
basic algorithm to provide for an accurate correction in the presence of strongly-absorbing 
aerosols, and addressing several MODIS-specific issues. 

L Validate/Improve the Basic MODIS, Algorithm using SeaWiFS 

We will begin the algorithm validation and improvement process using SeaWiFS imagery. 

As we have been able to successfully process SeaWiFS imagery (reformatted to MODIS) 
using our MODIS code, we can be assured that the implementation is correct. Thus, we 
will concentrate on the science of the correction algorithm by using SeaWiFS to validate the 
water-leaving radiances and indicate any potential problems and potential improvements. 
For example, we have found using SeaWiFS data acquired during the "initialization" cruise 
in January and February of 1998, that using the 670-865 nm bands (the "6-8" algorithm) 
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resulted in a somewhat better atmospheric correction than using the 765-865 nm bands (the 
"7-8" algorithm). This has important implications for MODIS, because it uses an algorithm 
that is similar to the 7-8 algorithm in SeaWiFS. However, the SeaWiFS band at 765 nm 
overlaps the 0 2 "A" absorption band, while the corresponding MODIS band (748 nm) 
avoids it. Thus, we need to know whether the better performance of the 6-8 compared to 
7-8 algorithm of SeaWiFS is due to the 0 2 "A" band, or a more fundamental problem that 
could have an impact on MODIS. 


2. Implement and Test the "Spectral Matching Algorithm" 

We will continue our basic study of the "spectral matching" algorithm. Unlike the basic 
MODIS algorithm, the spectral matching algorithm has the potential for atmospheric 
correction in the presence of strongly-absorbing aerosols. 

The spectral matching algorithm [Gordon et al., "Remote sensing ocean color and aerosol 
properties: resolving the issue of aerosol absorption," Applied Optics, 36, 8670-8684 
(1997); Chomko and Gordon, "Atmospheric correction of ocean color imagery: Use of 
the Junge power-law aerosol size distribution with variable refractive index to handle 

aerosol absorption," Applied Optics , 37, 5560-5572 (1998)] has already been implemented 
in an image processing environment. We will test its efficacy and document its 
performance (compared to the basic MODIS algorithm) using SeaWiFS imagery in regions 
in which the aerosol properties may be highly variable on a day-to-day basis, but the water 
properties are reasonably stable, e.g., the Middle Atlantic Bight (MAB) or the Tropical 
Atlantic and Caribbean. Our focus will be as complete an assessment as possible, so a 
decision can be made concerning its replacement of, or inclusion in, the basic MODIS 
algorithm. At a minimum, it will be used to provide a flag in the MODIS algorithm that 
will signal the probable presence of absorbing aerosols, and indicate that the quality of the 
derived products cannot be assured. However, the long-term goal is that it replace the 
basic algorithm. 


3. Optical Properties^ of WindzBlown Dust 

We need quantitative estimates of the optical properties of wind-blown dust, e.g., from 
Africa, to provide a proper model to effect atmospheric correction using the standard 
MODIS algorithm. 


The optical properties of wind-blown dust are required to effect atmospheric correction 
using the standard MODIS algorithm, and possibly the spectral matching algorithm as well. 
We have a unique opportunity now to use SeaWiFS imagery that was acquired over the 
Caribbean this year's particularly intense dust periods of July and August, simultaneous 
with our surface measurements of the aerosol vertical profile (using LIDAR) and aerosol 
optical thickness (using a CIMEL sun/sky radiometer) made from St. Johns, VI, to develop 
a model for African dust. 


4. MODIS-Specific Issues 

We need to address the detection and removal of thin cirrus clouds, methods for efficiently 
including earth-curvature effects, out-of-band corrections, BRDF effects on the diffuse 
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transmittance, correct for polarization sensitivity ofMODIS, and the efficacy of 
atmospheric correction for removal of the aerosol effect from the measurement of the 
fluorescence line height. 

Of the issues listed under this heading, the first that we will pursue is the MODIS out-of- 
band corrections. As we now have the MODIS relative spectral response (RSR) functions 
from MOST, we can complete their incorporation into the algorithms following the 
procedures described by Gordon (1995) ["Remote sensing of ocean color: a methodology 
for dealing with broad spectral bands and significant out-of-band response", Applied 
Optics, 34, 8363-8374 (1995)]. 

The second is incorporating the SBRS/MCST polarization-sensitivity data into the 
atmospheric correction module. This will be effected as described in Gordon, et al., (1997) 
["Atmospheric Correction of Ocean Color Sensors: Analysis of the Effects of Residual 
Instrument Polarization Sensitivity," Applied Optics, 36, 6938-6948 (1997)]. 


The third issue to be examined is the BRDF effect. As described below (In-water Radiance 
Distribution), we propose to use our measurements of the BRDF to develop a model that 
can be applied to MODIS imagery. This model will be used to address the BRDF on the 
diffuse transmittance. 

The forth issue in importance is efficiently including earth-curvature effects in the MODIS 
algorithm. Following Ding and Gordon ["Atmospheric correction of ocean color sensors: 

Effects of earth curvature," Applied Optics, 33, 7096-7106 (1994)] this will most likely be 
a modification of the look-up tables for the top-of-the-atmosphere contribution from 
Rayleigh scattering. However, before actually embarking on an implementation, we will 
examine SeaWiFS imagery at high latitudes to assess the impact of neglecting earth- 
curvature in the algorithm. 

Finally, because of uncertainty in the performance of MODIS Band 26, and because we 
need to assess whether our radiative transfer codes are sufficiently accurate to study 
removal of the aerosol effect from the measurement of the fluorescence line height, we will 
examine this issue and thin cirrus clouds only if time permits. 

An updated Algorithm Theoretical Basis Document (ATBD) for the Normalized Water - 
leaving Radiance algorithm will be prepared and submitted by April 30. 1 999. 


Whitecap Correction Algorithm 

Objectives and Proposed Activities (FY '99) 

Our basic goal for the rest of the project is to maintain experience in operating and 
maintaining the instrumentation in preparation for the validation phase of the contract, and 
to complete the analysis of the whitecap data acquired thus far. 
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5 Maintaining, Measurement Capability, 

We need to maintain our ability to make whitecap measurements. 

The basic objectives of the experimental portion of this task has been realized (acquiring 
whitecap radiometric data at sea), experimental work is being suspended until the validation 
phase, except insofar as the radiometer is being operated at sea when it is sufficiently 
important to do so, e.g., it was operated during the SeaWiFS Initialization Cruise (MOCE- 
4). This requires personnel capable of both maintaining and operating the instrumentation. 

6. Reduction and Analysis of Existing Data 

We need to complete the analysis of the data already acquired. 

We have a significant amount of whitecap data that is yet to be reduced and analyzed. In 
addition, we need to reanalyze the Tropical Pacific whitecap data because of the 
surprisingly low reflectance increase due to whitecaps that we measured there. This is a 
unique data set, as it was acquired in the trade winds with moderately high winds (8-12 
m/s) and practically unlimited fetch and duration. We have now developed an alternative 
method of analysis that we believe is more robust and will provide greater confidence in the 
results. Earlier this year, the SeaWiFS Project informed us that the present whitecap 
algorithm was causing the atmospheric correction algorithm to fail in the South Atlantic. 

We provided an algorithm adjustment (based on our earlier analysis of the Tropical Pacific 
whitecap data), and were informed by them that the algorithm appeared to be working 
much better after the adjustment. This example underscores the importance of a detailed 
examination of as much SeaWiFS imagery as possible prior to the launch ofMODIS. 

In-water Radiance Distribution (BRDF) 


Objectives and Proposed Activities (FY '99) 

During FY '99 our objectives for this task are: maintaining experience in operating and 
maintaining the instrumentation in preparation for the validation phase of the contract; 
acquiring more field data; and using the field data to develop a model of the BRDF as a 
function of the solar zenith angle and the water's chlorophyll concentration. 

7. Maintaining Measurement Capability 

We need to maintain our ability to make BRDF measurements. 

This requires personnel capable of both maintaining and operating the instrumentation, as 
well as reducing the data. 

8. Acquisition of More Field Data 

To build and thoroughly test a BRDF model, we need to acquire data over a wider range of 
chlorophyll concentrations. 

We will operate the Radiance Distribution Camera System (RADS) whenever the 
opportunity to acquire data over a wider range of chlorophyll concentrations presents itself. 
For example, we plan to participate in INDOEX in the Indian Ocean in January-February 
1999. This will also provide an independent data set for validation of the MODES 
algorithm using SeaWiFS. We will, of course, participate in the MODIS initialization 
cruise. 
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9. Build a BRDF Model 

We need to develop a model relating the subsurface radiance distribution (BRDF) to the 
chlorophyll concentration of the water 

The data acquired during the SeaWiFS initialization cruise showed that, at a solar zenith 

angle of ~37° (the mean encountered during the cruise) and a chlorophyll concentration of ~ 
0.1 mg/m 3 (the approximate mean encountered during the cruise), an error of ~5% in the 
normalized water-leaving radiance [L iv (A)] a , would be made if the measurement was carried 

out at nadir rather than at the angle appropriate to the viewing direction of the sensor. Since 
most investigators are only capable of measuring the in-water upwelling radiance at nadir, 
direct comparison of their measurements with MODIS (or SeaWiFS) data will result in an 
error at, or above, the level of error that is acceptable for the MODIS product. This error 
can be reduced by using a model (based on the chlorophyll concentration and solar zenith 
angle) to either correct the nadir measurement to the appropriate viewing angle, or to correct 
the MODIS product to provide the normalized water-leaving radiance at nadir. In either 
case, a model is required and we propose to build one that includes both elastic and 
inelastic (Raman) scattering, as it is evident that a considerable portion of the radiance in the 
blue-green (>10%) results from inelastic processes. In addition we will compare our 
measurements with the / / Q model of Morel [Morel and Gentili, "Diffuse reflectance of 
oceanic waters. III. Implication of bidirectionality for the remote sensing problem," 

Applied Optics, 35, 4850—4862 (1996)] (which is being used in other models in the 
community). Construction of such a model may require that the in- water light field be 
simulated with a radiative transfer code that includes polarization. 


Pre-launch/Post-Iaunch Atmospheric Correction Validation 

Objectives and Proposed Activities (FY '99) 

The original objectives of this task were fourfold: (1) study aerosol optical properties over 
the oceans to assess the applicability of the aerosol models proposed for use in the 
atmospheric correction algorithm; (2) measure the aerosol optical properties (including their 
vertical distribution) from a ship during the initialization/calibration/validation cruises; (3) 
determine how accurately the radiance at the top of the atmosphere can be estimated based 
on measurements of sky radiance and aerosol optical thickness at the sea surface (i.e., 
vicarious calibration); and (4) utilize data from other sensors that have achieved orbit 
(OCTS, POLDER, SeaWiFS ...) to validate and fine-tune the correction algorithm. We 
have obtained a significant amount of data toward (1), designed, constructed, or 
purchased, instrumentation to acquire data for (2), completed (3), and (4) was discussed 
under the first activity above (Atmospheric Correction Algorithm Development). 

Objectives for the post-launch validation phase are: maintain experience in operating and 
maintaining the instrumentation in preparation for the validation phase of the contract; 
complete analysis of data already acquired; and participate in the validation phase. 

10. Maintaining Measurement Capability 

We need to maintain our ability to make atmospheric measurements at sea. 
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This requires personnel capable of both maintaining and operating the instrumentation, as 
well as reducing the data. We will continue to keep the CIMEL operation in the Dry 
Tortugas, including the monthly maintenance checks. We plan to use the MPL, aureole 
camera, and sky camera during INDOEX in January-February of 1999, and during the 
MODIS initialization cruise. 

11. Complete Analysis of Existing Data 

We need to complete the analysis of the data acquired previously. 

We have been operating the CIMEL instrument in the Dry Tortugas almost continuously for 
several years. It worked well, largely due to the diligence of E. J. Welton in maintaining 
the instrument and the site. However, we have only scratched the surface in the data 
analysis. We have only looked at a few specific time periods, and need to examine the 
entire record. In particular, there are several periods when African dust was known to be 
present, and we can study its properties using the CIMEL data. Also, we have extracted 
specific days of the data set, believed to be dust, marine aerosol, or nonseasalt sulfate 

aerosols. On these days we are running our inversion method [Wang and Gordon, 
"Retrieval of the Columnar Aerosol Phase Function and Single Scattering Albedo from 
Sky Radiance over the Ocean: Simulations," Applied Optics, 32, 4598-4609 (1993)] for 

recovering the phase function, to compare with that used by Nakajima et al. ["Retrieval of 
the Optical Properties of Aerosols from Aureole and Extinction Data," Applied Optics, 22, 
2951—2959 (1983)], used in the Aeronet Network). This work will be continuing. 

Aureole and sky camera data acquired during the July Hawaii criuse are being analyzed, 
specifically for several locations while the cruise went near the volcanic plume to look at the 
retrieved size distribution of particulates in the plume. The aureole and sky cameras 
(including the sky polarization) were also operated during the MOCE-4 cruise, and these 
data are in the process of being reduced and analyzed as part of the SeaWiFS initialization. 

12. Post-launch Validation 


We will participate in MODIS post-launch validation. 

We will use the MPL, aureole camera, and sky camera during INDOEX in January- 
February of 1999. In this case, SeaWiFS will be a surrogate for MODIS in the algorithm 
validation program. We will also participate in the MODIS initialization cruise; however, 
its scheduling will be dependent on the MODIS launch. 

Detached Coccolith Algorithm 

Objectives and Proposed Activities (FY '99) 

The objectives of our FY 99 effort are to complete analysis of cruise work done to 
date, continue new field data in the Gulf of Maine and Gulf of Mexico, use these data to 
improve algorithm performance and to apply and validate the coccolithophore algorithm 
using SeaWiFS data. 

13. Processing of completed pre-launch cndses 

We have been collecting data on coccolith concentrations in several pre-launch 
cruises per year for the past several years. Most analyses have now been performed, and 
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post-cruise calibrations of the scattering measurements have been done (calibration 
standards are always run before each cruise, but periodic distilled water checks are always 
run, and must be checked against published values to verify calibration). While basic 
interpretation of the prelaunch data have been done, much more data interpretation is in 
order. Thus, we propose to devote half of our time in FY 99 to further reduction of the 
pre-launch data. This also will involve using SeaWiFS data of water-leaving radiance at 
550 nm to validate the algorithm, and comparing against shipboard measurements of 
suspended calcite. There was a coccolithophore bloom last year in the Bering Sea, in 
which the coccolithophore algorithm was implemented. Due to the fact that the bloom was 
extremely thick, we discovered that the look-up table must be re-run in order to apply the 
algorithm. This is now being done. 

14. New field efforts 

The other half of our effort will be devoted to pre- and post-launch cruises. We 
have several days of pre-launch cruises in the Gulf of Maine in September and October 
which will also have to be processed in the coming year. In April, 1999, we have 15 days 
of shiptime in the Gulf of Mexico between Tampa and Progresso, Mexico. We are 
including water-leaving radiance in our suite of continuous measurements, as sampled from 
a bridge-mounted Satlantic radiance sensor. We also will have 20 days of ship time next 
year in the Gulf of Maine in which we will collect coccolith concentrations, suspended 
calcite concentrations, inherent and apparent optical measurements. 

15. Improving algojnthm performance 

Now that SeaWiFS imagery can be run with MODIS software, we intend to 
compare the suspended calcite concentrations estimated from space-based measurements, 
with ship-board measurements. Specifically, we will begin with satellite-derived 
measurements from blooms, and work into nonbloom waters, in order to estimate the limits 
of sensitivity of the algorithm. Given that next year's Gulf of Maine cruises will be "post- 
launch", hopefully, we can test the algorithm without having to work through SeaWiFS. 
However, initial validation of the algorithms will be effected with SeaWiFS data. 


An updated Algorithm Theoretical Basis Document (ATBD)for the Detached Coccolith 
Concentration algorithm will be prepared and submitted by April 30. 1999. 


MODIS Ocean Discipline Group Validation Plan 

We will provide input to a revision of the Ocean Discipline Group's validation plan by 
March 30. 1999 . The plan will be in both narrative and view graph forms. 
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In this study, a new atmospheric correction algorithm capable of simultaneously retrieving 
aerosol and ocean optical parameters under the presence of both highly- and weakly- 
absorbing aerosols has been developed. It is suggested that the radiative properties of 
realistic aerosols can be well simulated with those resulting from the Junge power-law 
aerosol models. The use of the latter makes it possible to vary the atmospheric radiative 
properties continuously through a variation of the aerosol parameters. The atmosphere 
is assumed to consist of two, plane parallel and horizontally homogeneous, layers with 
the Fresnel reflecting bottom boundary. The radiative properties of the ocean water are 
assumed to be those of Case 1 waters. A system of non-linear equations is constructed 
for the radiances detected by a multi-band remote sensor and, is subsequently solved using 
non-linear optimization procedures. The algorithm’s performance has been studied with 
simulated test data. It is shown that the aerosol single scattering albedo (n7 0 ) and the 
pigment concentration (C) can be excellently retrieved to within 6% and 10% respectively 
even under the presence of the instrument calibration errors. However, because of significant 
differences in the scattering phase functions for the test and power-law distributions, large 
error is possible in the estimate of the aerosol optical thickness. The positive results for 
the pigment concentration C suggest that the detailed shape of the aerosol scattering phase 
function is not needed for the atmospheric correction of ocean color sensors. The relevant 
parameters are the single scattering albedo and the relative spectral variation of the optical 
depth. The vertical distribution of aerosols and a spectral variation of the aerosol’s refractive 
index have adverse effects on the accuracy of retrievals. Fortunately, these cases are easily 
identifiable in the course of non-linear optimization procedure as long as the data “fit” 
objective function SlsQ becomes relatively large, i.e., 5-6%, in contrast to less than 1% 
in other test cases. The algorithm was incorporated into the SeaWiFS image processing 
system SeaDAS. The results demonstrate that the algorithm’s performance is superior to 
the NASA Standard Atmospheric Correction algorithm. A significant advantage of the new 
approach is that realistic multicomponent aerosol models are not required for the retrieval 
of C. 
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I. INTRODUCTION 


The subjects of the atmospheric particulate turbidity and ocean color have been re- 
ceiving considerable attention recently. An impetus for the renewed interest is scientific: 
the role of aerosol and phytoplankton in the global carbon cycle and possibly climate. 1 

Marine phytoplankton are the microscopic single-celled ocean plants, and axe at the 
base of the marine food chain. Phytoplankton contain chemical chlorophyll-a and other 
pigments that absorb light, thus, providing the energy needed for photosynthesis. The rate 
of photosynthesis in phytoplankton is regarded as primary productivity. The study of the 
primary productivity and the global distribution of phytoplankton can provide answers to 
its role in the global climate forcing and carbon cycle. 2 ’ 3 Phytoplankton pigments absorb 
light mainly in the red (670 nm) and in the blue (440 nm) range of the spectrum, and 
reflect the light in the green (550 nm). By measuring the light exiting the ocean in the blue 
and green (the “ocean color”), the phytoplankton pigment concentration cam be determined 
(pure water reflects the light in this region of the spectrum). 

The passive remote satellite sensors provide an efficient amd cost-effective way to study 
the variation of marine phytoplamkton on a global scale. The light signail at the top of the 
atmosphere is comprised of two parts resulting from the radiative processes within: 1) the 
atmosphere itself (90-100%), and 2) the ocean (0-10%). Thus, the atmospheric radiative 
effects need to be removed. This problem is regarded as atmospheric correction. 

Atmospheric correction is complicated by many factors. Chief among them are aerosols 
(a suspension of fine liquid and/or solid particles in a gas). The aerosols possess not only 
a complicated chemical character but also diverse physical and geographical properties: 
shape, size distribution, time variability, location, etc. 4 ’ 5 There is no way to predict their 
properties a priori. 


1 
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The proof-of-concept ocean-color mission was the Coastal Zone Color Scanner (CZCS) 
on Nimbus-7 (1978-1986). 6,7 It demonstrated that the concentration of marine phytoplank- 
ton could be estimated from measurements of the radiance leaving the top of the atmo- 
sphere. This led to the development of more sophisticated ocean color sensors such as the 
Sea- Viewing Wide-Field-of-View color sensor (SeaWiFS), 8 the moderate resolution imag- 
ing spectrometer MODIS, 9 the polarization and directionality of the earth’s reflectances 
instrument (POLDER), 10 etc.). 

The first atmospheric correction algorithm that was capable of retrieving water-leaving 
radiances 6 ’ 7 was developed for the Coastal Zone Color Scanner, and utilized a simple single- 
scattering algorithm. Gordon and Wang 11 substantially improved it later including multiple- 
scattering effects and was successfully applied for SeaWiFS. This algorithm provided water- 
leaving radiances with the required accuracy for typical marine atmospheres. It failed, 
however, in cases when the atmosphere contained highly absorbing aerosols. 

Recently, Gordon et al. 12 proposed an algorithm for Case 1 waters 13 (i.e., waters whose 
optical properties are controlled by phytoplankton and their immediate detrital material) 
that worked resonably well for both weakly and highly absorbing aerosol and was able to 
obtain similtaneously the pigment concentration. It systematically varied the ocean pigment 
concentration and aerosol models, selected a pre-set number of best fits (in an r.m.s. sense) 
and computed the average of the model parameters taken to be the results of retrievals. The 
candidate models were derived from the bimodal log-normal size distributions and refractive 
indices provided by Shettle and Fenn 14 as a function of relative humidity. The water-leaving 
radiance was computed on the basis of the semianalytic radiance model of ocean color of 
Gordon et ai ls 

For the CZCS imagery correction Zhao and Nakajima 16 also proposed an algorithm 
that was able to determine aerosol and ocean properties similtaneously. They employed 
models with Junge power-law size distributions with a single index of refraction and used 
water-leaving ratios instead of water-leaving radiance itself. They showed good retrievals 



3 


with simulated test data created by use of bimodal aerosol size distributions with the same 
index of refraction and concluded that the power-law size distribution was an adequate 
approximation to the now favored bimodal aerosol distributions. 

This dissertation extends the work of Gordon et al. 12 through utilizing the Junge power- 
law models with variable refractive index and including non-linear optimization procedures 
to determine similtaneously aerosol and ocean water properties. We begin in Chapter 2 with 
the basic description of radiative transfer through the atmosphere-ocean coupled system. 
Chapter 3 provides the details of the proposed atmospheric correction algorithm. The 
performance of the newly developed algorithm is tested with simulated data is provided in 
Chapter 4. This algorithm has been incorporated into the SeaDAS image processing system 
for SeaWiFS and the results are demonstrated in Chapter 5. Chapter 6 concludes this work. 
The quasi-Newton minimi zation theory concepts are outlined in the Appendix. 



II. RADIATIVE TRANSFER IN ATMOSPHERE-OCEAN SYSTEM 


In the scalar approximation (the polarization ignored), the propagation of light through 
the plane-parallel, horizontally homogeneous atmosphere is governed by the Radiative 
Transfer Equation (RTE) 17,18 : 

C ° S 9 dr = ~ L ( r] ft + J ( r; 6 1 ft' (2.1) 

where 9 and <f> are the photon’s polar and azimuth angles, the optical thickness r = 
Jq c{z*)dz where c(z') is the extinction coefficient of the atmosphere at the geometrical 
depth z (assuming that z = 0 at the top of the atmosphere), and L{t\ 9, <f>) is the radiance, 
/(r; 9 , <j>) is the source function defined as 

J(r]d,<f>) = ~ r do' f 2 * d<f>' sinfl' L(t-, O' , <j>') , (2.2) 

Jo Jo 

where P(r; 9 1 , 0, <f>) is the single scattering phase function from the direction {9\(f> t ) to 

(0,<f)). The single scattering albedo tzr 0 = &(t)/c(t), where 6(r) is the scattering coefficient 
of the atmosphere at the optical depth r. In order to obtain a unique solution in the 
atmosphere two boundary conditions are needed: (1) at the top of the atmosphere 

L{ 0; 0, <f>) = E 0 *(cos 9 — cos 9 0 )S(<f> — <£ 0 )> 0 < ^r/2, (2.3) 

where Fo is the extraterrestrial solar irradiance, 9 0 and <j > o are the sun’s zenith and azimuth 
angles respectively, and (2) at the botom of the atmosphere 

L(n-,e r ,4> r ) = J r(9 it 4>i -> o r ,<f> r ) dsi it ( 2 . 4 ) 

where subscripts “r” and “i” stand for “reflected” and “incident” respectively, and r(9i,<f>i 
— > 9 ri <f> r ) is the surface’s bidirectional reflectance distribution function (BRDF). 

To represent the atmosphere-ocean coupled system we adopt a widely used plane-parallel 
horizontally homogeneous two-layer model for the atmosphere with the Fresnel reflecting 


4 
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bottom boundary where the top layer contains a gas of molecules (for simplicity, referred 
to as a Rayleigh layer ; ’ ) and the second bottom layer contains only aerosol particles unless 
specified otherwise ( u aerosol layer”). Radiance entering the ocean is backscattered into the 
atmosphere in accordance with the semi- analytic model of Gordon et. al. 15 


Let us introduce a dimensionless quantity named reflectance 


p(X) = 


ttL(A) 

i'o(A) cos 6 0 ’ 


(2.5) 


where L is the radiance, Fo is the extraterrestrial solar irradiance, 0 O is the solar zenith 
angle, A is the wavelength. The radiance at the top of the atmosphere (TO A) can be written 
then as (Figure 2.1) 


Pt Pr pa d“ Pra Pg *i” t>p w -j- tp WC} (2.6) 

where p t is the total reflectance measured by a remote sensor, p r is the reflectance resulting 
from the multiple scattering by the molecular (Rayleigh) gas only (i.e., without being 
scattered by aerosol), p a is the reflectance resulting from the multiple scattering by the 
aerosol only (i.e., without being scattered by Rayleigh gas), p ra is the interaction term that 
represents the reflectance resulting from scattering by both Rayleigh gas and aerosol, p g the 
reflectance resulting from the sun glitter, p wc is the reflectance resulting from whitecaps, 
p w is the water-leaving reflectance and, finally, t is the diffuse transmittance function of the 
atmosphere. 

As long as the ocean color sensors are equipped with tilting mechanisms, the term p g 
can be omitted. Assuming further that p wc = 0 the Eq. (2.6) simplifies to 


Pt = pr + pa + Pra + tPun (2*7) 

Our objective (and that of the atmospheric correction in general) is to obtain p w {^) from 
the values of p t (\). Below, a focus will be placed upon the computational aspects of each 
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of the terms shown in the Eq. (2.7). 


ATMOSPHERIC CORRECTION 

PROBLEM STATEMENT 



Figure 2.1. Atmospheric correction model. 
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Computation of 
Assuming that 

- the air molecules are considered as molecules of a single gas; 

- the scattering coefficient per unit mass of the air molecules /c ai > and the acceleration 
of gravity g do not vary with the height; 

- the molecular weight fraction of various atmospheric air molecules does not vary with 
the height, 

a simple expression for the optical thickness r r may be obtained as the function of the 
optical thickness r r0 corresponding to the standard surface pressure P 0 = 1013.25 mb and 
the ambient pressure P 

T r = T r 0 J^, (2.8) 

Hansen and Travis further derive an empirical relationship for r r o at STP: 

T r0 = 0.008569 A~ 4 (1 + 0.0113 A~ 2 + 0.00013 A~ 4 ), (2.9) 

where A is the wavelength in [/mi]. It should be noted, however, that more accurate 
computations may be performed without making the above assumptions, 19 but the deviation 
of the optical thicknesses of the Rayleigh gas for the atmospheric conditions provided by in 
the “TJ.S. Standard Atmosphere Supplements, 1976” 20,21 never exceeded 1% comp ar atively 
to those obtained with the equations (2.8)-(2.9) . Figure 2.2 demonstrates the results of the 
similar computations. 

In the absence of the polarization of the incident radiation and the molecular anisotropy, 
the angular distribution of the light scattered by the molecular gas is given by the Rayleigh 
phase function P(0) 


P(Q) = 3/4 (l + cos 2 (0)) 


( 2 . 10 ) 
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where © is the scattering angle. 


For a given sun- viewing geometry, 9 0 ,9 and the ambient pressure P, Gordon et al. 22 
show that the Rayleigh reflectance p r (A) in the visible and near-infrared can be computed 
through 


Pr(A; 6, 6 0 , &<f>) 


1 — exp(— r r (A)/ cos 9) 
1 — exp(— r r o(A)/ cos 9) 


p r o(X; 9, 9 0 , A <j>), 


( 2 . 11 ) 


where p r o(A; 9, 9q, A <f>) is the reflectance of the Rayleigh gas at STP, and A <j> is the viewing 
azimuth angle relative to the sun. 



Figure 2.2 Rayleigh optical thickness r r (A) for the U.S. 
Standard 1976 and four Supplements 1976 atmospheres 


The Eqs. (2.8) - (2.11) along with the Radiative Transfer equation constitute a complete 
set of equations sufficient to determine the radiative properties and the reflectance p r (A) of 
the Rayleigh gas. It is worth pointing out that the impact of the Rayleigh gas reflectance 
p r (A) on the resulting total reflectance at the TOA, p t ( A), decreases rapidly from the blue 
to the near infrared (NIR) spectral range, more specifically, p r (NIR) comprises just about 
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30-40% of p t (NTR) depending mostly on the sun-viewing geometry, in contrast to 70-80% 
in the blue. 

Computation of pA A) 4- p-^fA) 

The aerosol is composed of particles of extremely complex physical and chemical properties 
that depend on many factors like the aerosol origin, location, the time of a day, the 
wind speed, humidity, etc. They come in different shapes, forms and sizes. Their total 
unpredictability makes the problem of the atmospheric correction nearly untreatable. 

We start with a simplification that the particulates are uniform spheres which is a valid 
assumption for the pure water aerosols found in abundance over large bodies of water, like, 
oceans. The theory of the scattering of electro-magnetic plane waves from dielectric spheres 
was originally developed by Mie 23 in 1908. More modem treatments of the theory along 
with further developments are given by van de Hulst 18 and Kerker. 24 The theory for the 
electro-magnetic scattering from spheroidal particles is provided by Mishchenko. 25 

In this work, we will adopt the Mie theory in the development of the atmospheric 
correction algorithm. There exists an extensive list of literature on the subject of the 
Mie theory, hence, below we provide only those properties that are directly used in the 
development of our algorithm. 

The aerosol is characterized by its index of refraction 

m(A) = m r (A) — zmi(A), (2*12) 

where m r (A) and mi(A) are respectively the real and imaginary part of the index of 
refraction. D ’Almeida, Koepke and Shettle 4 provide the refractive indexes for the most 
common aerosols. They axe partially reproduced in the form of a diagram in the Fig. 2.3. 

It can be observed that the imaginary part m;(A) of the index of refraction of water- 
soluble and dust aerosols have only a slight dependence on the wavelength, while the rest 
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are strongly dependent on the wavelength. Note, however, that neither oceanic (or sea-salt) 
nor sulfate (mostly existing over land) aerosol is ever found in its “pure” form and come 
as a mixture of various aerosols (for example, the maritime aerosol consists of 99%, by 
number, of the water-soluble and 1% of the oceanic components as defined by Shettle and 
Fenn 14 ). In the atmospheric correction, we will make the assumption that the refractive 



Wavelength, nm 

Figure 2.3. Imaginary part of the index of refrac- 
tion mi (A) for selected aerosols: water-soluble, dust, 
oceanic, maritime-mineral, sulfate. Based on the data 
provided by d’Almeida, Koepke and Shettle, 1991 

index is independent of the wavelength. 

Another aspect that determines the optical properties of aerosols is the spectrum of 
particle sizes. The smelliest particles (typically called Aitken nuclei ) range in size from 
10“ 3 /xm to 10” 1 /xm and have a little effect on the optical phenomena. The condensation 
nuclei or large particles are important, on the other hand. They range in size from 10 1 /xm 
to 1 fan. Giant particles are those greater then 1 ^m in radius (see, Elterman 26 ’ 2 for the 
detailed classification). Deirmendjan (1960)“ 8 parameterized the aerosol size distribution 
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with the following equation 

n ( D ) = = aD a e ~ bD1 , (2.13) 

where n(D) is the size frequency at the diameter D, dN(D ) is the number of particles with 
diameters between D and D + dD and a, a, 6 , and 7 are positive constants determined 
by the type of the aerosol. Besides the fact that the parameters in (2.13) had a readily 
interpretable physical meaning, they fit fairly well with measurements on natural water 
clouds and aerosols. However, we will adopt another parameterization that proves to be 
convenient when dealing with aerosols that are a mixture of different types of aerosols, or 
modes : the log-normal size distribution 

M 

n P) = E ni ( i) )’ ( 2 - 14 ) 

t=i 


where n*(D) is the size frequency at the diameter D for the Vs component (mode) of the 
size distribution and M is the total number of modes. For the individual component n t *(Z?) 


we have 


log.(10)v^2> 


exp - 


logio {D/Di) 


where 07 is the standard deviation, Ni is the total number density of the ith component, 
Di is the mode’s diameter. On the basis of the Eq. (2.15) Shettle and Fenn 14 and later 
d’Almeida, Koepke and Shettle 4 developed an extensive number of aerosol models having 
optical properties encountered under realistic conditions (some of them will be considered 
in later chapters in more detail). Having derived the scattering a sca and absorption <r a bs 
cross-sections from the Mie theory for individual particles of the diameter D, we can further 
compute quantities that directly govern the radiative transfer such as the scattering and 
absorption coefficients, k aca and k abs , respectively (note that the total extinction coefficients 
is defined to be the sum k cxt = k sca + k abs = c), e.g., 


ksca= \ 2 cr sca (D)n{D)dD J (2.16) 

J £>1 

where we explicitly restricted the range of diameters from D\ to D 2 in the light of the 
above discussion. The coefficients k abs and k cxt are defined by the similar expressions. The 
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single-scattering albedo (representing the fraction of the radiance removed from a beam of 
light by scattering) follows immediately from fc-coefficients 


GJq — ksca/kexti 


(2.17) 


Figures 2.4 and 2.5 provide the single-scattering albedo t«7o(A) and the extinction coefficient 
k e x *(A) as the function of wavelength A. These diagrams reveal a remarkable features of 
functional behaviours of ur 0 (A) and k ext ( A). They are very smooth functions of the argument 
A, and in particular, k ext ( A) demonstrates apparent log-linear relationship over the spectral 
range 400 — 900 nm. 



Wavelength, nm 

Figure 2.4 Single-scattering albedo tu'o(A) as the function 
of wavelength A [ft m], for selected aerosols: water-soluble, 
dust, oceanic, maritime- mineral, sulfate. Based on the 
data provided by d’Almeida, Koepke and Shettle, 1991 


The angular distribution of the scattered light is called the phase function P(0), where 
© is the scattering angle. It is also computed from Mie theory in a manner similar to 
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k S ca ^ the Eq.(2.16). Examples of P(0) are provided in the next chapter. These optical 



400 500 600 700 800 900 

Wavelength, nm 

Figure 2.5 Extinction coefficient k ext { A) as the function 
of wavelength A [/xm], for selected aerosols: water-soluble, 
dust, oceanic, maritime-mineral, sulfate. Based on the 
data provided by d’Almeida, Koepke and Shettle, 1991 

parameters, k exty U7 0 , and P(0) provide sufficient information to solve the radiative transfer 
equation (RTE) within the atmosphere for p a + p ra . Once computations of p a + Pra are 
carried out for discrete values of 0 , and r a , it is necessary to make some approxima- 

tions in order to interpolate between computed values in an image processing environment. 
In our work we extensively use the approach outlined by Gordon and Wang. 11 In particular, 
they have shown that the sum of the terms p a (A) + p ra { A) at a given sun- viewing geometry 
is a smooth function of the aerosol optical thickness r a (Figure 2.6 demonstrates such a 
relationship for four of Shettle and Fenn’s Maritime (M80), Tropospheric (T80), Coastal 
(C80), and Urban (U80) aerosol models, all at the relative humidity 80%) and can be 
adequately described by the four coefficients a, 6, c, d as 

Pa{ A) + Pra(A) = a(A) T a ( A) + 6(A) t J ( A) + c( A) rl( A) + d( A) r^( A), (2.18) 


Those coefficients are stored in the form of look-up tables (LUT). We will return to this 
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issue in later chapters where a further storage reduction technique is described in more 
detail. We conclude this subsection with a note that the dependence of p 0 (A) + p ra (A) on 



Aerosol optical thickness, t a 

Figure 2.6 Reflectance p a (^) + /Va(A) as a function 
of the aerosol optical thickness r a for the solar zenith 
angle 0 O = 20° and viewing angle 0 — 1°. Based on 
the RTE computations 


r r ( A) through the r r (A) dependence on pressure is extremely small and is neglected. 29 


Computation of the diffuse transmittance t( A) 


The water-leaving radiance at the top of the atmosphere (TOA) is related to 

the water-leaving radiance L w (£ v ) at the bottom of the atmosphere through the diffuse 
transmittance £(£„) function defined as 

t(L) = ( 2 - 19 ) 

where is a unit vector directed from the sea-surface to a sensor and £*£,(&,) is the water- 
leaving radiance just above the sea- surface. While £(&,) is a weak function of the aerosol 
radiative properties, its impact, generally speaking, cannot be neglected. It also has been 
shown 30 that t(£ v ) is a function of the angular water-leaving radiance distribution. In this 
work, two assumptions have been examined: 
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1. The aerosol influence on £(f v ) is ignored and the water-leaving radiance just beneath 
the sea surface is uniform, and 

2. The aerosol influence on t(£ v ) is included. 


Assumption 1 leads to a simple relation for t(£ v ) empirically derived by Gordon and Clark. 31 

+f\ \ ^ r ( IriA) , _ \ ( 1 , 1 M 


t(X) = exp — 


1 1 
cos 9 0 cos 9 V 


where r r (A) is the Rayleigh optical thickness, 9q and 9 V axe the solar zenith angle and the 
sensor viewing angle with respect to the zenith, r 0z { A) is the optical thickness of the ozone 
layer located just above the Rayleigh layer on TOA. 


Computations of the diffuse transmittance t(£ v ) under the Assumption 2 are more 
involved. Yang and Gordon (1998) replace (2.20) by 

t{L) = A(e v ) (2.21) 

where A(9 V ) and B(9 V ) are coefficients resulting from a fit to directly computed values of 
t(f„). Note that these coefficients are the functions of the sensor polar single only. Further, if 
the aerosol is non-absorbing, t(| v ) is only a weak function of the aerosol optical thickness r a , 
and independent of the aerosol vertical structure, and if the aerosol is absorbing t(£ v ) is still 
independent of the aerosol vertical structure. In this work, Eq.(2.21) is used. This means 
we include no coupling of the ocean and atmosphere through t(A), and no bidirectional 
effects of L w (£ v ) on t( A). 


Computation of the water-leaving reflectance 


Gordon et a/. 15,3t defined the normalized water-leaving reflectance as 

V(l-/9)(l-p)iZ " 
m 2 Q( 1 — rR) 


( 2 . 22 ) 


where p is the Fresnel reflectance of the sea surface for normal incidence; p is the Fresnel 
reflection albedo of the sea surface for irradiance from the Sun and sky; R is the irradiance 
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reflectance just beneath the sea surface (i2 = E u /Ed, where E u end Ed are the upwelling 
and downwelling irradiances just beneath the sea surface respectively); m is the index of 
refraction of water; and Q is the ratio of the upwelling radiance (at nadir) to the upwelling 
irradiance toward the zenith ( Q equals tv for a totally diffuse radiance distribution); r is the 
water-air interface reflectance for totally diffuse irradiance. 

We will consider briefly the ocean optical properties for Case 1 waters only, 32 " 34 
i.e., waters whose optical properties are controlled by phytoplankton and their immediate 



Wavelength, nm 

Figure 2.7 Normalized water-leaving reflectance Z^A) as 
the function of wavelength A [/mi] for two pigment pigment 
concentrations 0.05 and 1.0 mg/m 3 , and three values of the 
scattering parameter 6°: 0.12, 0.30, 0.45 m" 1 . 

detrital material. According to the bio-optical model developed by Gordon et al., 15 p w can 
be computed from two parameters: 6° relating the scattering coefficient to the pigment 
concentration, and the pigment concentration itself (pigment concentration is defined as 
the stun of the concentrations of chlorophyll a and phaeophytin a). The value of 6° was 
empirically determined 34 to be in the range from 0.12 to 0.45 m" 1 having the average value 
of 0.30 mT 1 for the pigment concentration C ranging from 0.05 mg/m 3 to 10 mg/m 3 for 
Case 1 waters. The Figure 2.7 provides sample computations of the normalized water- 
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leaving reflectance for pigment concentrations 0.05 mg/m 3 , 1.0 mg/m 3 and three values of 
the scattering parameter 6°. Note the strong dependence of [p w ]w on C. 



HI. SPECTRUM OPTIMIZATION ATMOSPHERIC CORRECTION ALGORITHM 


Of all developed aerosol models for spherically shaped particulates the Junge power-law 
aerosol models are the simplest, and are described by the size distribution 

K , Do < D < D x 


dN 


7d = \ k (*/*) 


v+1 


(3.1) 


D\ < D < D 2 

0 , D < Z?o, jD > Z?2 

where diV is the number of particles per unit volume with diameters between D and D+dD, 

K is the normalization constant such that J 0 °° dN = 1. It is sometimes more convenient 
to use the so called volume size distribution dV(D) that can be obtained from the particle 
size distribution through the transformation of variables 


dV(]nD) = -D z dN(D), 

6 

Shown in the Fig. 3.1 are the volume size distributions for v = 2.0, 3.0, 4.0 


(3.2) 



Particle Diameter. p.m 

Figure 3.1. Volume size distribution for the Junge power- 
law aerosols with v — 2.0, 3.0, 4.0. 
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An outstanding feature of Junge power-law size distributions is that all their optical 
properties are controlled by just one parameter, v. Van de Hulst 18 shows that if the particle 
refractive index is independent of A, D x = 0 and D 2 = oo, then k ext oc A~ 2+l/ and, vj and 
P(@) are independent of A. For the values of D 0l Di and D 2 used here, these facts are 
still approximately true. Figure 3.2 demonstrates the fraction of the extinction coefficients 
c(A)/c(A = 865 nm) as the function of wavelength (c = k ext ). Curves corresponding to 
different parameters of v are easily identifiable and, in fact, are almost independent of 
the refractive index m(m r ,m t ). Independent of the refractive index is also the amount of 
forward scattering as can be seen from Fig. 3.3. The amount of backward scattering is a 
strong function of the imaginary part mi of the index of refraction m. As expected, the 
shape of the phase function P(0) doesn’t vary much with wavelength (up to 30% at 120°). 



400 520 640 760 880 

Wavelength, ntyi 

Figure 3.2. The fraction c(A)/c(865) as the function of 
wavelength for the Junge power-law aerosols with v = 

2. 0(0. 5)4.0. 





v = 4.0 
v = 3.0 
v = 2.0 


rT^r**0:040r* 


Figure 3.3. The scattering phase function P(0) as the 
function of the scattering angle © for the Junge power-law 
aerosols with v = 2.0, 3.0, 4.0 and selected m, at A = 865 


Figure 3.4. The single scattering albedo z&q as the func- 
tion of the imaginary part of the index of refraction m t - for 
the Junge power-law aerosols with selected parameters of 
A, v and m r . 









21 


Far more complex (see Fig. 3.3) is the behaviour of the single scattering albedo zu 0 
on its parameters. The strongest variation appears to be with the imaginary part of the 
index of refraction in fact, it is an exponentially decreasing function of m,- and a weak 
function of the real part m T of m. At higher m* the single scattering albedo becomes 
a strong function of wavelength. The dependence on v is moderately strong and highly 
non-linear. 

In the Chapter 2 we have made a note of a smooth, almost linear, behaviour of the 
reflection function p a (A) + Pra(^) on the aerosol optical thickness r a (see Eq. (2.15)). The 
purpose of the thesis is to exploit this property utilizing Junge power-law size distributions 
described above. We begin with the following 

Proposition. A reflectance field /> t (A) — /? r (A) produced by essentially 
any aerosol and any ocean Case 1 waters can be adequately regenerated 
by the reflectance field resulting from the Junge power-law model aerosols 
described by a set of parameters and the ocean bio-optical 

model described by a set of parameters {C, 6°}. Furthermore, this set of 
values {i/, m r , mi, r a , C, 6°} is unique. 

If the above conjecture is true, and the parameters are obtained one way or another, we 
end up having a set of useful aerosol and ocean water characteristics. We now study this 
proposition using simulations. 

Let us introduce some simplifying notation: 

V, m r , 77li, T a ) = p a (AjJ I/, m r , mi, T a ) -f- Pra(^j j ^~a) j 

paw(Xj', C,b°) = v,Tn r ,Tm,T a ) p w {\y, C,b°), (3.3) 

pw( A;; C,b°) = t(Aj) Pw ( A i; C,b°), 

where t M (Xj ; t/, m r , m ( , r a ) and t(\j) are the diffuse transmittance functions in the sense of 
equations (2.21) and (2.20) respectively, p w {Xj; C,b° ) is the water-leaving reflectance. 
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Typically, ocean color sensors measure the radiance in a finite number N of pre-selected 
spectral bands (the SeaWiFS color sensor, for instance, is supplied with N=8 bands in the 
visible and NTR spectral ranges). In general, N non-linear equations in the six variables 
{u, m r ,mi,T a ,C, b 0 } can be written for N bands (see, Eq. (2.4): 

Pa(*i] v,m r ,Tni,T a ) + p aW ( A t ; C,b°) = p^) - p^), 

Pa( A 2 ; v, m r , mi, r a ) -(- p a \v{\ 2 -, C,b°) = p t {\ 2 ) - Pr(A 2 ), 

(3-4) 

/>4(Ajv; i',m r ,m i ,r a ) + p a w(Ajv; C,b°) = p t (\ N ) - p r (\ N ). 

Note that each of variables {v,m r ,mi, C, 6°)} is subject to constraints: 

0 < v (min) < v < v (max) , 

0 < m;(min) < m; < max) , 

0 < m r (min) < m r < m r ( max) , (3.5) 

0 < C (min) < C < C (max) , 

0 < 6°(min) < 6° < 6° (max) . 

Here we assume C'(min) = 0.1 mg/m 3 , C(max) = 1.5 mg/m 3 , 6° (min) = 0.12 m -1 , £>°(max) 
= 0.45 m -1 , and i/(min) = 2.0, i/(max) = 4.5, m r (min) = 1.33, m r (max) = 1.50, mi(min) 

=0.0, m t (max) = 0.04. The ranges of 6° and v are consistent with oceanic and atmospheric 

measurements. 


Now, we choose a least-square objective function (LSQ) 5 L sq(^, m r , rrii, r a , C, b°) such 

that 


Slsq^, m r ,mi,T a ,C,b ) — 


N 


N 


i - 


, 77i{ , T a ) Hh Pa\v(^ji C/fb ) 

Pt(^j) ~ Pr{^j) 


(3.6) 


where the numerator within the square brackets is to be calculated and the denominator 
is provided by a sensor (assuming that the reflectance p r (Aj) can be calculated exactly for 
any wavelength given the ambient surface pressure, as has been discussed previously in the 
Chapter 2). The LSQ objective function serves as a measure of “goodness” in the fit of 
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Pa + PaW to p t - p r . In other words, at the solution point {i/*, ra*, m* , r*, C*, 6°*} 5 L sq 
should attain its minimal value subject to constraints in (3.5). 

A number of versatile methods exist to find such a solution, 35 although they might 
differ appreciably as to the convergence rate, the memory requirements, and the ability to 
coverge to the solution. In particular, methods belonging to the so-called quasi-Newton 
class 36 (see also the Appendix) that take advantage of the information contained in the 
Hessian matrix of a function (i.e., the square matrix of the second partial derivatives of 
the function calculated at a point) in calculating the step size and the direction of search 
appear to be adequate for our purposes. Powell 37 proves the convergence to a minimum for 
the Davidon-Fletcher-Powell (DFP) 35,35>33 variable metric class method (which also holds 
for Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm) on a function with properties 
outlined above in Eqs.(3.4)-(3.5). However, to the best of our knowledge, the quasi-Newton 
class methods do not guarantee, in general, a convergence to a minimum. In this study 
we used a slightly modified DZXMWD routine from the commonly available International 
Mathematical & Statistical Library 39 (IMSL), which is based on the Harwell Library routine 
VA10A and incorporates a quasi-Newton method. To assure the convergence to a minimum 
in our non-linear optimization procedure a set of starting search points is taken, some of 
which are discarded after performing a few iterations based on their relative values. It 
has been observed that the resulting solution was the same (within a preset tolerance) 
whatever the number of starting search points chosen, meaning that the convergence to a 
globed minim um was probably achieved, so it is believed as sufficient to use just one starting 
point in solving the system (3.4), and the result is the solution we want. 

In the meantime, the system (3.4) subject to constraints (3.5) is solved in two different 
ways: 1) with a reduction of the number of variables, and 2) directly solving (3.4). Both 
ways have advantages and disadvantages on their own. Though this will be discussed later, 
here, we note that the system with a reduced number of variables becomes more stable, it 
achieves the desired accuracy of computations faster, and affirmatively provides the global 
minimum of the solution. 
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Solution 1. 

It can be noted that the water-leaving radiances in the NIR bands (i.e., at A = 765 
nm, 865 nm for SeaWiFS) are generally extremely small comparatively to the radiances 
reflected by the atmosphere (indeed, they constitute, typically, no more then 0.1% of the 
total radiance measured by a sensor at the TOA). As a consequence, we disregard terms 
P wipin', C,6°), such that p a w(A 7 ; C,b°) = p a w(A 8 ; C, 6°) = 0 in the system (3.4). 


v - Approximation 



v 


Figure 3.5. Relationship between £'(765,865) and v 
(solid curve) averaged over all combinations of four values 
of r a (865) (0.1, 0.2, 0.3, 0.4), six values of m* (0.0, 0.001, 
0.003, 0.010, 0.040), and two values of m r (1.334 and 1.50), 
with 0o = 20° and 9 = 45.9°. Dashed curves represent 
the upper and lower envelopes of the relationship, and the 
darker region represents the standard deviation of the var- 
ious cases about the mean for the given value of u. 


Gordon 40 has shown that the spectral variation of p a (^NiR) + Pra(^ N ir) in the NIR 
depends mostly on the aerosol size distribution and only slightly on the refractive index, 
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thus the NIR bands can be used to determine the size distribution parameter v. In fact, 
introducing a measured value 




Pa(Aj) ~f~ Pra(Aj) 
Pa{Xl) + P ra ( A*)’ 


(3.7) 


where X a and A / are the shorter and the longer of the two NIR bands available for the 
sensor, and averaging it over all model values of {r a , m r , ra*}, 




Nm r + Nm i + N Tc 


iV mr ,iV miI iV ra 

mV, m[ r \ r^), (3.8) 


t,r,fc=l 


where JV mr is the total number of preselected values m for the real part of the refractive 
index, N mr is the total number of preselected values for the imaginary part of the 
refractive index, and N Ta is the totcil number of preselected values for the optical 
thickness r a . We can form an empirical relationship between 1 1 and the size distribution 
parameter v. It is shown in Figure 3.5. 


Once e r (\ si \i) is computed from the sattelite imagery, it is a straightforward matter 
to determine the Junge power-law parameter i/*, where the * indicates the estimated value. 
From the tabulated values of p a (765) + p ra (765) and /> a (865) + />ra(865) for each model from 
the set {m r ,m t *}, we derive the aerosol optical thickness r a at 865 nm (Fig. 2.5) through 
the Newton- Raphson iteration procedure 41 . This iteration proved to be fast and accurate. 
Note that the opticad thicknesses r a (A) for other wavelengths can be found (see, Fig. 3.2) 
because 

_c(A)_ _ Tq(A) (3 .9) 

c(865) r 0 (865) ’ v 

where c(A) is the extinction coefficient computed through Mie theory for each set {z/, m r ,mj} 

and properly interpolated for the calculated value of v * . 

In building the system (3.4) for the bands not containing the NIR wavelengths we 
have to assume some function for the diffuse transmittance of the water-leaving reflectance 
through the atmosphere. This has been discussed previously in the Chapter 2. Both 
Assumptions 1 and 2 (Eqs. (2.20) and (2.21) respectively) are easily introduced in these 
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computations. However, we will choose the simpler one of the two mentioned (i.e., we will 
use theEq. (2.20). This will mean, in fact, that the reflectances [/>a(A)+p ra (A)] and [pw] are 
completely separable and unrelated to each other. While, in general, such an assumption 
is not justified, Yang and Gordon 30 showed that the diffuse transmittance t(X) for small 
optical thicknesses is only a weak function of the aerosol optical properties. 

Thus, the total number of variables {i/, r a , C , 6°} is reduced by two {i/, r a }. The 

system (3.4) no longer contains equations corresponding to the NIR bands and is re-written 
below for clarity: 

Pa(Ai; i/,m r ,mi,T a ) + p w ( A x ; C, 6°) = p t (A x ) - p r {\i), 

Pa{ A 2 ; v,m r ,Tni,T a ) + p w (\ 2 ] C, b°) - p t (X 2 ) - />r(A 2 ), 

(3.10) 

Pjt(Ajv_jv N i R ; v,Tn r ,mi,T a ) + pw{ Ajv; C,b°) = p t (\ N ^N mR ) — p r (X N ~ NmR ). 

Now, applying the quasi-Newton iteration subject to constraints (3.5) with the objective 
function 5 L sq from (3.6), we obtain the required solution (7*, 6° }. A 

detailed graphical flow-chart in the Figure 3.6 demonstrates all the steps performed in 
this solution. 

Solution 2. 


In this case, the quasi-Newton iteration is applied to the full system of equations (3.4) 
(including the NIR bands) subject to same constraints (3.5) and the objective function S LS q 
from (3.6) 

Discussion. The ultimate test of any correction algorithm is its stability and robustness. 
This was the main reason for introducing two solutions above (though in a sense similar). 
It is surely possible under realistic conditions to encounter aerosols with optical properties 
considerably different from those described by the simple Junge power-law aerosol models. 



SPECTRUM: OPTIMIZATION ALGORITHM: Flow-Chart 


1. Determine V 
X= 765, 8 65 


u 

«'* 

a 

sr (765,865) 

R. (765) 
R- (865) 

I “ 



O.* 

l 


2. Determine x a co 
X = 865 nm 


this curve is obtained 
from pre-computed LUT's 
for all X, m £ , m r 
and given viewing geometry 


£'(765,865) 

GIVEN 


GIVEN 

curw It obtained 


Wtw/n^ gmommtry 


c / c(865) 


3. Determine x a at 
X = 412, 443, 490, 510, 555, 670 

c / c(865) = t/t(865) 

these values are computed from 
model* s scattering parameters 
for all X, m m r 


4. Determine p a + p ra 

X = 412, 443, 490, 510, 555, 670 

r/zzj’ curve is obtained 
from pre-computed LUT’s 
for all X, m £ , m r 
and given viewing geometry 


5. Solve system of NLE with. LSQ-objective 


this curve is interpolated. 



GIVEN 



guessydi^ie 


solution 


Note: the procedure might be thought as follows: 
"guess" a value for C, bO arid calculate p w 
then find (mi,mr) such that P a ■+■ P ra + Pw = P t “ 


Figure 3.6, Atmospheric Correction Flow-Chart 
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Yet, we require a good performance with a little noise, with the retrieved parameters as 
close to the real ones as possible. 

An apparent drawback of the Solution 1 is its inability to obtain the exact parameters 
for pseudo data generated using the Junge power-law for the aerosols. The reason is the 
estimate of ^(A*, A/) in (3.8). Since the deviation of e cxa ct from? might reach up to 3%, it is 
not possible to expect better accuracy in obtained values for {z/, m r , m*}. On the other hand, 
the usage of (3.8) makes the algorithm’s performance extremely stable in situations when 
the test aerosol is different from Junge power-law aerosol models. The retrieved parameters 
are obtained with the standard deviations never exceeding 15-20%. It also provides for a 
better global optimal solution. As will be shown in the later chapters, the algorithm is 
extremely robust when dealing with the experimental errors such as the sensor calibration 
errors, when dealing with the aerosol vertical structures, etc. The Solution 2, on the other 
hand, computes the parameters {i/, m r ,mj} “exactly” but its results are less stable, the 
standard deviation reaching sometimes 30%. 

While both the solutions are extremely time-efficient (about 2 sec for an image pixel 
on AlphaServer 2100), the Solution 1 is preferable to the Solution 2. However, the Solution 
2, because of its direct approach, is more flexible, and can be easily extended to include, 
for example, the aerosol vertical structures. Here, the Solution 2 is used mostly to test the 
correctness of the computer code used in the implementation of the Solution 1. 



IV. ALGORITHM’S PERFORMANCE WITH SIMULATED DATA 


In the following subsections we generate and test the algorithm with simulated values 
for the radiances measured by the pseudo-sensors at the TOA. They will provide: 

- a theoretical background on test aerosols or probable errors to be expected with actual 
color sensors; 

- the results for the retrieved atmospheric-ocean parameters; 

- a brief discussion of results. 

A. Pseudo pA A) for Junge power-law test aerosol models 

The purpose of this test is to examine the correctness of all interpolation procedures 
as well as the performance of the optimization itrerations. 

A test data-set was simulated by computing p t ( A) — p r (A) for the Junge power-law 
aerosol models with m = 1.50 — i0.002 (characterized by low absorption) and an optical 
thickness r a (865) = 0.2 at A = 865 nm. Parameters for the ocean water properties were 
chosen to be C = 0.1 mg/m 3 and 6° = 0.3 m~ l (the average value of 6° for Case 1 waters 15 ). 
Calculations were performed for seven sun- viewing geometries with the solar zenith angles 
9 0 = 0°, 20°, 40°, 60° and, the sensor viewing at 9 V = 1°,45° with respect to nadir and 
the azimuth angle <f> v = 90° (omitting the glitter at 9 0 = 0 °,9 V = 1°). The atmospheric 
correction scheme of the Solution 2 (namely, the direct solution of the system (3.4)) was 
employed in this case. The retrived results for parameters {i/*, ra*, zUq, t*(865), C *, 6° } are 
shown in the Table 4.1. 
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Table 4.1 


a) 


\ 

0 

® 0 

0° 

20° 

O 

60° 

, 0 

V* 


2.07 

2.04 

2.05 

1 

m * 


1.500 - i 0.003 

1.500 - i 0.002 

1.500 -i 0.003 

45° 

V* 

2.06 

2.06 

2.05 

2.04 

m * 

1.500 - i 0.003 

1.453 - i 0.002 

1.500 - i 0.003 

1.500 - i 0.002 


b) 


\ 

e 

® o 

0° 

20° 

o 

O 

60° 

, o 

E9I 


0.910 

0.930 

0.929 

1 



0.24 

0.20 

0.20 

45° 


0.924 

0.932 

0.929 

0.934 


0.20 

0.22 

0.20 

0.20 


c) 


X 

0 

\ 9 ° 

0° 

O 

O 

<N 

o 

© 

© 

o 

o 

* 0 

C * 


0.101 

0.099 

0.098 

1 

b o 


0.30 

0.30 

0.30 

43° 

C * 

0.098 

0.098 

0.098 

0.098 

bo 

0.30 

0.30 

0.30 

0.30 


Table 4.1 Retrieved aerosol parameters for the Junge power-law 
test model with m r = 1.50, m; = 0.002, v — 2.0,r a (865) = 0.2 
and, the ocean water pigment concentration and the scattering 
parameter C = 0.1 mg/m 3 and b° = 0.3 m -1 respectively, for 
seven sun- viewing geometries. The correct value of ctq is 0.944. 
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Figure 4.1. a) Retrieved the total, aerosol and water- 
leaving reflectances for the Junge power-law test model 
with m r = 1.50, m; = 0.002, v = 2.0,r a (865) = 0.2 and, 
the ocean water pigment concentration and the scattering 
parameter C = 0.1 mg/m 3 and b° = 0.3 m -1 respectively, 
for the sun-viewing geometry 9o = 20°, 0 — 1°. b) respec- 
tive retrieval errors 
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Figure 4.2. a) Retrieved the toted, aerosol and water- 
leaving reflectances for the Junge power-law test model 
with m r = 1.50, mi = 0.002,1/ = 2.0,r a (865) = 0.2 and, 
the ocean water pigment concentration and the scattering 
parameter C = 0.1 mg/m 3 and 6° = 0.3 m” 1 respectively, 
for the sun-viewing geometry 0 o = 60 °, 9 = 1°. b) respec- 
tive retrieval errors 
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The retrieved and given reflectances along with the error deviations are plotted in 
Fig. 4.1 and 4.2 for two sun-viewing geometries, 0 O = 20°, 9 = 1° and 0 Q = 60°, 0 = 1°, 
<t>o = 90°. Apparently, the aerosol and ocean water parameters are retrieved excellently: the 
maximum deviation for the single scattering albedo zu 0 never exceeded 3%, the deviation for 
the pigment concentration C never exceeded 2%. As expected in the previous chapters, the 
retrieval of the aerosol optical thickness r a ( A) is the most unstable (with respect to both the 
sun- viewing geometries and the other parameters), though stable enough to be acceptable. 
The results demonstrate that the adopted interpolation and optimization procedures work 
satisfactorily, so that other test cases should be considered. This also demonstrates that 
the algorithm was properly coded. 


B. Pseudo p*( A) for M80. C80. T80, U80 test aerosol models 

The nomenclature M80, C80, T80, U80 stands for the Maritime, Coastal, Tropospheric 
and Urban models of Shettle and Fenn, 14 all at the relative humidity of 80%. The dis- 
tinguished feature of these models is that they closely approximate optical properties of 
realistic aerosols and are characterized by the bimodal size distribution (see Eq. (2.12) with 
i = 2). 


Table 4.2: Characteristics of the M80,C80,T80,U80 test aerosol models. 



Size Distribution 

Refractive Index 

Aerosol Model 

Ni 

Di 

< Ti 

412 nm 

865 nm 

M80 

0.990000 

0.06548 

0.35 

1.446— i 3.309e— 3 

1.436— i 6.107e— 3 


0.010000 

0.63600 

0.40 

1.359 — i 5.165e— 9 

1.348— i 1.381e— 6 

C80 

0.995000 

0.06548 

0.35 

1.446— t 3.309e— 3 

1.436— i 6.107e— 3 


0.005000 

0.63600 

0.40 

1.359— i 5.165e— 9 

1.348— i 1.381e— 6 

T80 

0.990000 

0.06548 

0.35 

1.446— i 3.309e-3 

1.436— i 6.107e— 3 

U80 

0.999875 

0.07028 

0.35 

1.423— i 3.473e— 2 

1.414— i 3.412e— 2 


0.000125 

1.16200 

0.40 

1.415— t 3.151e— 2 

1.406-i 3.095e— 2 


The size distribution parameters and values of the refractive index at two wavelengths 
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A = 412 and 865 nm are tabulated in the Table 4.2. Figure 4.3 compares the volume size 
distribution dV/dlnD for the test models M80, C80, T80, U80 and the Junge power- 
law models with v = 2.0, 3.0, 4.0. All volume size distributions are normalized such 
that Io° = 1 cm’ 3 . Figure 4.4 demonstrates the spectral variation of the aerosol 
optical thickness with respect to that at A = 865 nm. We note a strong similarity of 
this almost log-linear behaviour with that of the Junge power-law models (see Fig. 3.2). 


Table 4.3: The single-scattering albedo zjq of the M80,C80,T80,U80 test aerosol models. 


Wavelength, nm 

M80 

C80 

T80 

U80 

412 

0.9924 

0.9884 

0.9758 

0.7823 

443 

0.9929 

0.9890 

0.9761 

0.7825 

490 

0.9935 

0.9898 

0.9763 

0.7820 

510 

0.9938 

0.9901 

0.9764 

0.7810 

555 

0.9936 

0.9897 

0.9737 

0.7805 

670 

0.9944 

0.9906 

0.9715 

0.7738 

765 

0.9939 

0.9896 

0.9634 

0.7621 

865 

0.9934 

0.9884 

0.9528 

0.7481 


From Table 4.3 it can be observed that the spectral variation of the single- scattering 
albedo H7 0 (A) is unsubstantial, being the highest for Tropospheric T80 models, again, a 
property which is similar to the Junge power-law models. Among all the models the Urban 
at the relative humidity of 80% possesses the lowest value of the single-scattering albedo 
(tsrj 80 (865) = 0.748 comparatively to the M80 zj™ 80 (865) = 0.993). The U80 model is 
referred to as the highly absorbing aerosol in contrast to M80, C80, T80 models which are 
usually referred to as the weakly absorbing aerosols. 

The test-data set [pt(A) — p r (A)] was created for M80, C80, T80, U80 aerosol models 
having the optical thicknesses r a (865) = 0.1, 0.2, 0.3. The ocean water optical properties 
were chosen to be 6° = 0.3 m” 1 and three values of the pigment concentration C = 
0.1, 0.5, 1.0 mg/m 3 . 
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Figure 4.3. Volume size distributions for the M80, T80, 
U80 test aerosol models and the Junge power-law test 
model with v = 2.0, 3.0, 4.0. 


Figure 4.4. Spectral variation of the optical thickness 
r a (A)/r 0 (865) for the M80, C80, T80, U80 test aerosol 
models. 
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The solar zenith angles 9 0 were 0°, 20°, 40°, 60° and the viewing angles 9 = 1°,45° with 
respect to nadir (total of the seven sun- viewing geometries since 9 0 = O°,0 = 1° case was 
excluded as long as it would be strongly influenced by the Sun-glitter). 

The retrievals were performed using the approach outlined in the previous chapter as 
Solution 1, where SeaWiFS bands 7 and 8 (at A = 765,865 nm respectively) were used to 
determine the Junge power-law parameter v and the optical thicknesses r a (A). Five bands 
1 through 5 (at A = 412,443,490,510,555 nm) were used in the non-linear optimization 
procedure, Eqs. (3.4)-(3.6). 

The detailed analysis will be carried out for the aerosol models M80, T80, U80. The 
C80 aerosol test model behaves similarly to the M80 model. It should be noted that the 
M80 aerosol phase function is highly peaked in the forward direction because of the high 
presence of the giant particles in its particle size distribution. The U80 aerosol is highly 
absorbing, and the T80 aerosol model is characterized by a high reflectivity which is the 
consequence of the dominant presence of particles in the region of 0.5 /xm characterized 
by the high scattering efficiency. As was mentioned before both M80 and T80 models are 
weakly- absorbing aerosols. 

Figures 4.5 through 4.13 have been brought in to demonstrate an unessential depen- 
dence of the retrieved results on the sun- viewing geometry. The tests are shown for the 
case with r a (865) = 0.2, the pigment concentration C = 0.5 mg/m 3 , the water scattering 
parameter 6° = 0.3 m" 1 , and three sun- viewing geometries 9$ = 20°, 9 = l° 7 <f> = 90°; 
9 0 = 60°, 9 = 1°,<£ = 90°; 0 O = 40°, 9 = 45°, <j> = 90°. In the latter geometry the optical 
path through the aerosol layer is the longest among all the geometries considered. Figures 
4.7, 4.10, 4.13 clearly indicate that for this geometry the retrieved aerosol-ocean parameters 
sure slightly poorer then for the other geometries. In all cases the aerosol reflectance are 
best retrieved for bands 7 and 8 (the ones that were used to obtain the Junge power-law 
parameters u and r a (A)) and the worst near the band 6 (this band was not used during the 
solution of non-linear equations). Outstanding is the fact that the algorithm provides the 
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results equally well for the highly absorbing U80 and the weakly absorbing M80 aerosols. 
The average value of the objective function 5 LS q with respect to the sun- viewing geometries 
was 0.24% for M80, 0.17% for U80, and 0.31% for T80 aerosol models. These numbers were 
typical while working with all the above test models and indicate that the fit to the total 
aerosol reflectance p t ( A) - p r { A) was excellent. 

It has been noticed in the previous chapter that the normalized water-leaving re- 
flectance is the highest for waters with the low pigment concentration, i.e., C = 0.1 mg/m" 3 
(see Fig. 2.6). As a consequence, these waters’ parameters should be the easiest to determine 
compared to waters with the high pigment concentration content, e.g., C = 1.0 mg/m" 3 . 
This is validated for two aerosol optical thicknesses r a (865) = 0.1, 0.3: r a (865) = 0.1 
in Figures 4.14 - 4.15 for M80, Figures 4.16 - 4.17 for T80, and Figures 4.18 - 4.19 for 
U80; and r a (865) = 0.3 in Figures 4.20 - 4.21 for M80, Figures 4.22 - 4.23 for T80, and 
Figures 4.24 - 4.25 for U80. The results are provided for just one sun-viewing geometry 
0q = 60°, 9 = 1°, (f> = 90°. In U80 case, for example, for the optical thickness r a (865) = 0.1 
the retrieval of the pigment concentration C gave the error 0% for C = 0.1 mg/m -3 and 
2.6% for C = 1.0 mg/m" 3 . 

The same figures demonstrate that with the increase of the aerosol optical thickness 
the aerosol parameters are retrieved better (the reflectance of aerosol increases), while the 
ocean parameters are retrieved more poorly. The ocean parameters are better retrieved 
for low r a . For example, for M80 model (Fig. 4.14 and 4.20) and the ocean pigment 
concentration C = 0.1 mg/m -3 the retrieved of C degraded from 4% for r a (865) = 0.1 to 
8% for r a (865) = 0.3 while the single-scattering albedo tzr 0 (865) improved from about 6% 
to 2% and r a (865) from about 30% to 7%. 

The results for the single-scattering albedo ru*o(865), the ocean pigment concentration 
C, and the water scattering parameter 6° are summarized in Tables 4.4 - 4.6. All values 
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were averaged over the sun-viewing geometries as 


A = 


Ni Nj 


w+^-)fer,tr 


EEMV’ 1 ) 


(4.1) 


where A stands for one of the parameters tcr 0 (865), (7, b° and, 9q , 0 are the sun and viewing 
zenith angle. 


From Table 4.4 we observe that the error of the retrieval of the ocean pigment concen- 
tration C was on the average 3.6% for M80, 4.5% for C80, 10.3% for T80, and 5.1% for U80 
test aerosol models. The algorithm’s handling the highly absorbing aerosol case (U80) is 
outstanding. The single-scattering albedo t^ 0 is remarkably well retrieved also (Table 4.5). 
For the weakly absorbing aerosols M80, C80 the average error was 2.7% and, for the highly 
absorbing aerosol U80 it was 2.6%. Notable is the stability of the algorithm’s performance. 
The objective function S LS q was on the average 0.27% for M80, 0.28% for C80, 0.36% for 
T80, and 0.15% for U80 aerosol test models with its typical average standard deviation of 
20% over all cases examined. 

An intrinsic drawback of the algorithm is its inability to retrieve reasonably the aerosol 
optical thickness r a . The correct single-scattering phase function is a primary condition 
for that. The “true” and “retrieved” phase functions for the test cases M80 and U80 are 
shown in Figures 4.26 and 4.27 respectively. Meanwhile, Tables 4.4 and 4.6 suggest that 
the retrieval of the ocean water parameters C and 6° is only a weak function of the single 
scattering phase function. 

We end up this subsection with a note that only the first five out of six bands in visible 
were used in constructing the systen of non-linear equations (3.4), (3.6) subject to con- 
straints (3.5). The addition of the Band 6 (670 nm) did not improve the result of retrievals. 
Indeed, the variation of the water-leaving reflectance p^(C, 6°) at this wavelength is small 
(but not negligible), thus downgrading the overall time-performance of the algorithm. 
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Figure 4.5 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-80 test model having 
the optical thickness r a (865) = 0.2 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.5 mg/m 3 and 6° = 0.3 m” 1 respectively, for the sun- 
viewing geometry 0 O = 20°, 9 = 1°,0 = 90°. b) respective 
retrieval errors 
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Figure 4.6 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-80 test model having 
the optical thickness r a (865) = 0.2 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.5 mg/m 3 and 6° = 0.3 m” 1 respectively, for the sun- 
viewing geometry 6q = 60°, 9 = l°,<f> = 90°. b) respective 
retrieval errors 
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Figure 4.7 a) Retrieved the total, aerosol and water- 
leaving reflect ances for the Maritime- 80 test model having 
the optical thickness r a (865) = 0.2 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.5 mg/m 3 and b° = 0.3 m” 1 respectively, for the sun- 
viewing geometry 9q = 40°, 9 = 45°, <f> = 90°. b) respective 
retrieved errors 
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Figure 4.8 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Tropospheric-80 test model 
having the optical thickness r a (865) = 0.2 and, the ocean 
water pigment concentration and the scattering parameter 
C = 0.5 mg/m 3 and 6° = 0.3 m _1 respectively, for the sun- 
viewing geometry 9 0 = 20°, 9 = 1°,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.9 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Tropospheric-80 test model 
having the optical thickness r a (865) = 0.2 and, the ocean 
water pigment concentration and the scattering parameter 
C = 0.5 mg/m 3 and b° = 0.3 m _1 respectively, for the sun- 
viewing geometry Oq = 60°, 9 = 1°,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.10 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Tropospheric- 80 test model 
having the optical thickness r a (865) = 0.2 and, the ocean 
water pigment concentration and the scattering parameter 
C - 0.5 mg/m 3 and 6° = 0.3 m” 1 respectively, for the sun- 
viewing geometry 9 0 = 40°, 9 = 45°, <j> = 90°. b) respective 
retrieval errors 
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Figure 4.11 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.2 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.5 mg/m 3 and 6° = 0.3 m" 1 respectively, for the sun- 
viewing geometry 9 0 = 20°, 9 = 1 °,0 = 90°. b) respective 
retrieval errors 
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Figure 4.12 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.2 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.5 mg/m 3 and b° = 0.3 m -1 respectively, for the sun- 
viewing geometry 9q = 60°, 9 = 1 °,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.13 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.2 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.5 mg/m 3 and 6° = 0.3 m~ l respectively, for the sun- 
viewing geometry 9q = 40°, 9 = 45°, (j> = 90°. b) respective 
retrieval errors 
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Figure 4.14 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-80 test model having 
the optical thickness r a (865) = 0.1 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.1 mg/m 3 and 6° = 0.3 m~ l respectively, for the sun- 
viewing geometry 0 Q = 60°, 0 = 1 °,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.15 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-80 test model having 
the optical thickness r a (865) = 0.1 and, the ocean water 
pigment concentration and the scattering parameter C = 
1.0 mg/m 3 and 6° = 0.3 m^ 1 respectively, for the sun- 
viewing geometry 9q = 60°, 0 = 1°,0 = 90°. b) respective 
retrieval errors 
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Figure 4.16 a) Retrieved the totcil, aerosol and water- 
leaving reflectances for the Tropospheric-80 test model 
having the optical thickness r a (865) = 0.1 and, the ocean 
water pigment concentration and the scattering parameter 
C = 0.1 mg/m 3 and 6° = 0.3 m” 1 respectively, for the sun- 
viewing geometry 9 0 = 60°, 0 = 1°,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.17 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Tropospheric-80 test model 
having the optical thickness r a (865) = 0.1 and, the ocean 
water pigment concentration and the scattering parameter 
C = 1.0 mg /m 3 and 6° = 0.3 m~ l respectively, for the sun- 
viewing geometry 9q = 60°, 9 = l°,<j> = 90°. b) respective 
retrieval errors 
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Figure 4.18 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.1 and, the ocean water 
pigment concentration and the scattering parameter C — 

0.1 mg/m 3 and 6° = 0.3 m~ l respectively, for the sun- 
viewing geometry 9q = 60°, 9 = 1°,0 = 90°. b) respective 
retrieved errors 
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Figure 4.19 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.1 and, the ocean water 
pigment concentration and the scattering parameter C = 
1.0 mg/m 3 and 6° = 0.3 m" 1 respectively, for the sun- 
viewing geometry 9q = 60°, 6 = 1°,<£ = 90°. b) respective 
retrieval errors 





60 .G? 0 =1.02? <j> = 90.0? #bands = 5, b 


p f - p r 



X 0 ° Retrieved Reflectances 


,-rj a " a -A. 

Pra + Pa 


0 . 020 - * 

\ 

tPw X 


C = 0.1 (mg/m 3 ) C ’ = 0.092 (mg/n? ) 


^*=0.3 
CO 0 = 0.993 
b°= 0.30 


^; = 0.32 
C0;= 0.970 
b° = 0.29 

S d = 0.41% 


Wavelength, nm 


total reflectance 
aerosol reflectance 
water-leaving reflectance 


Wavelength, nm 

Figure 4.20 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-80 test model having 
the optical thickness r a (865) = 0.3 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.1 mg/m 3 and b° = 0.3 m” 1 respectively, for the sun- 
viewing geometry &o = 60°, 0 = 1°,0 = 90°. b) respective 
retrieval errors 
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Figure 4.21 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-80 test model having 
the optical thickness r a (865) = 0.3 and, the ocean water 
pigment concentration and the scattering parameter C = 
1.0 mg/m 3 and b° = 0.3 m"* 1 respectively, for the sun- 
viewing geometry Oq = 60°, 9 = 1°,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.22 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Tropospheric- 80 test model 
having the optical thickness r a (865) = 0.3 and, the ocean 
water pigment concentration and the scattering parameter 
C = 0.1 mg/m 3 and 6° = 0.3 m -1 respectively, for the sun- 
viewing geometry 0q = 60°, 9 = 1°,<£ = 90°. b) respective 
retrieval errors 





57 




400 450 500 550 600 650 700 750 800 850 900 


Wavelength, nm 



Wavelength, nm 

Figure 4.23 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Tropospheric-80 test model 
having the optical thickness r a (865) = 0.3 and, the ocean 
water pigment concentration and the scattering parameter 
C - 1.0 mg/m 3 and 6° = 0.3 m~ l respectively, for the sun- 
viewing geometry 0q = 60°, 9 — 1 °,<£ = 90°. b) respective 
retrieval errors 
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Figure 4.24 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.3 and, the ocean water 
pigment concentration and the scattering parameter C = 
0.1 mg/m 3 and 6° = 0.3 m’ 1 respectively, for the sun- 
viewing geometry 0q = 60°, 9 = 1°,0 = 90°. b) respective 
retrieval errors 
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Figure 4.25 a) Retrieved the total, aerosol and water- 
leaving reflectances for the Urban-80 test model having 
the optical thickness r a (865) = 0.3 and, the ocean water 
pigment concentration and the scattering parameter C = 
1.0 mg/m 3 and 6° = 0.3 m” 1 respectively, for the sun- 
viewing geometry 9q = 60°, 0 = 1°,<£ = 90°. b) respective 
retrieval errors 
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RETRIEVAL OF OCEAN’s PIGMENT CONCENTRATION 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (&I80, C80, T80, U80). 

Also , provided are the standard deviations over viewing geometries divided by the mean (S d ) 
as well as the deviations from given parameters ( D ). 
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6 . 1 % 

2 . 10 % 

x a (865) =0.20 

0.086 

1.9% 

14.0% 

0.425 

9.9% 


0.941 

7.0% 

5.90% 

-c a (865) =0.30 

0.090 

5.7% 

10.0% 

0.412 

6.8% 

17.6% 

0.874 

6.7% 

12 . 6 % 



C; mg/m 3 

0.100 

0.500 

1.000 


result 

s d 

D 

result 

s d 

D 

result 

El 

D 

U 80 

CD 0 = 0.748 

i 

■ 

x a (865) =0.10 

0.100 

0.8% 

0.00% 

0.507 



1.028 

1.9% 

2.80% 

t a (865) =0.20 

0.104 



0.526 



1.053 

3.2% 

5.30% 

x a (865) =0.30 

0.110 



0.552 

5.8% 

10.4% 

1.070 

4.0% 

7.00% 


Table 4.4 
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RETRIEVAL OF AEROSOL’s SINGLE SCATTERING ALBEDO 

Tabulated are mean values of retrieved C0 0 for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80 , C80, T80, U80). 

Also , provided are the standard deviations over viewing geometries divided by the mean (S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/m 3 

0.100 

0.500 

1.000 


result 

EM 

D 

result 

El 

D 

result 

s d 

D 

M80 

C0q = 0.993 

x a (865) =0.10 

0.938 

2.5% 

5.50% 

0.968 



0.972 

2.2% 

2.12% 

x a (865} =0.20 

0.954 

2.2% 

3.93% 

0.971 

2.4% 

2.22% 

0.974 

2.5% 

1.91% 

x a (865) =0.30 

0.963 

1.8% 

3.02% 

0.976 

1.9% 

1.22% 

0.979 

2.0% 

1.41% 



C; mg/m 3 

0.100 

0.500 

1.000 


result 

EM 

D 

result 


D 

result 

mm 

D 

C80 

co 0 =0.988 

x a (865) =0.10 

0.943 

2.3% 

4.56% 

0.970 

1.7% 

3.20% 

0.972 

1.7% 

1.62% 

x a (865) =0.20 

0.960 

2.1% 

2.83% 

0.977 

1.7% 

1.11% 

0.979 

1.8% 

0.91% 

x a (865) =0.30 

0.969 

1.9% 

1.92% 

0.981 

1.6% 

0.71% 

0.983 

1.5% 

0.51% 



C; mg/m 3 

0.100 

0.500 

1.000 


result 

is 

D 

result 

s d 

D 

result 

EM 

D 

T80 

cOq = 0.953 

x a (865) =0.10 

0.842 

4.9% 

11.6% 

0.867 

4.5% 

9.02% 

0.872 

4.4% 

8.50% 

x a (865) =0.20 

0.877 

1.9% 

8.00% 

0.891 

2.0% 

6.51% 

0.892 

1.9% 

6.40% 

x a (865) =0.30 

0.902 

1.6% 

5.35% 

0.914 

1.0% 

4.09% 

0.918 

1.0% 

3.67% 


!!!!!!!!®!!!!i 

C; mg/m 3 

0.100 

0.500 

1.000 


result 

Sd 

D 

result 

Sd 

D 

result 

s d 

D 

U80 

COq =0.748 

x a (865) =0.10 

0.744 

4.7% 

0.54% 

0.745 

4.3% 

1.40% 

0.745 

4.4% 

0.40% 

x a (865) =0.20 

0.770 

0.8% 

2.94% 

0.769 


2.81% 

0.770 



x (865) =0.30 

0.788 

1.7% 

5.35% 

0.777 

1.2% 

3.88% 

0.773 




Table 4.5 
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RETRIEVAL OF WATER SCATTERING PARAMETER 

Tabulated are mean values of retrieved b ° for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80 , T80, 1180). 

Also, provided are the standard deviations over viewing geometries divided by the mean (S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/m 3 

0.100 

0.500 

1.000 



result 

El 

D 

result 

El 

D 

result 

s d 

D 

M80 

x a (865) =0.10 

0.294 

1.1% 

2.00% 

0.274 

4.1% 

8.67% 

0.273 

4.4% 

3.20% 

COg = 0.993 

•t a (865) =0.20 

0.289 

1.7% 

3.67% 

0.258 

7.3% 

14.% 

0.255 

8.2% 

15.0% 


x a (865) =0.30 

0.286 

1.8% 

4.67% 

0.249 

8.2% 

17.0% 

0.246 

9.4% 

18.0% 



C: mg/m 3 

0.100 

0.500 

1.000 


result 

El 

D 

result 

s d 

D 

result 


D 

C80 

co 0 =0.988 

x a (865) =0.10 

0.293 

1.6% 

2.33% 

0.273 

4.9% 

3.20% 

0.271 

5.3% 

9.67% 

■c a (865) =0.20 

0.289 

2.2% 

3.67% 

0.255 

8.0% 

15.0% 

0.250 

9.0% 

16.7% 

t a (865) =0.30 

0.287 

2.0% 

4.33% 

0.247 

7.4% 

17.7% 

0.239 

8.7% 

20.3% 



C; mg/m 3 

0.100 

0.500 

1.000 


result 

s d 

D 

result 

El 

D 

result 

s d 

D 

T80 

C0 0 = 0.953 

x a (865) =0.10 

0.278 

6.5% 

7.33% 

0.237 

t— ‘ 

OO 

JP 

21.0% 

0.234 

19.% 

22.% 

■c a (865) =0.20 

0.270 

6.4% 

10.0% 

0.200 

23.% 

33.3% 

0.193 

23.7°/ 

b 36.% 

x a (865) =0.30 

0.288 

5.1% 

4.00% 

0.201 

20.% 

33.0% 

0.186 

25.% 

38.0% 



C; mg/m 3 

0.100 

0.500 

1.000 


result 

s d 

D 

result 

Sd 

D 

result 

mm 

D 

U80 

C0 0 = 0.748 


0.297 


1.00% 

0.298 

1.8% 

0.67% 

0.298 

1.7% 

0.67% 

t a (865) =0.20 

0.302 

0.5% 

0.67% 

0.316 



0.319 

2.4% 

6.33% 





0.350 



0.354 

8.1% 



Table 4.6 



















































































































































mm 


1.0 ( mg/m 3 ) 


v = 2.41 


Figure 4.26 Comparison of the true (solid curve) and 
retrieved (dashed curve) single-scattering phase functions 
at 865 nm for M80 test aerosol model 


1.0 (mg/m 3 ) 


v = 3.26 


Figure 4.27 Comparison of the true (solid curve) and 
retrieved (dashed curve) single-scattering phase functions 
at 865 nm for U80 test aerosol model 
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C. Pseudo fo r M80. C80. T8Q. U80 test aernsnl mnrUU 
with added sensor calibration error 

To assess the effect of the sensor calibration errors in p t { A) on the algorit hm perfor- 
mance we adopt the approach outlined by Gordon. 42 The pseudo-reflectance data p t ( A) 
measured by a remote sensor were simulated for models M80, C80, T80, U80, as in the 
previous subsection, and a pseudo-calibration error added to all bands as follows: 

p'tW = Pt(^) [1 + ct(A)], (4.2) 

where a(A) is the fractional error in p t ( A), p' t { A) is the pseudo-reflectance at the TOA 
measured by a sensor. 


Table 4.7: Values of the Residual Radiometric Uncertainty 
after Effecting an In-Orbit Calibration Adjustment 


A,, nm 

Uncertainty, % 

412 

0.3 

443 

0.5 

490 

0.8 

510 

1.0 

555 

1.5 

670 

2.0 

765 

3.0 

865 

5.0 


Gordon 42 has demonstrated an in-orbit calibration procedure for which the residual 
calibration error in all the visible bands will have the same sign as in NIR bands and, the 
magnitude of the error will progressively decrease from NIR to visible bands roughly as the 
Rayleigh reflectance p r increases (which is calculated exactly and has the greatest impact on 
the toted reflectance p t in the blue). Assuming ~ 5% calibration error at 865 nm, following 
Gordon 42 the values of the residual radiometric uncertainty a(A) for all bands after effecting 
the vicarious calibration procedure are provided in the Table 4.7. We consider both positive 
and negative calibration errors. 
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The retrieval of the ocean pigment concentration has been summarized in the Table 
4.8 (positive calibration errors) and Table 4.9 (negative calibration errors) where parameter 
values were averaged over all sun-viewing geometries (Eq. (4.1)). The “true” and retrieved 
reflectances are shown for two particular cases with added positive calibration error: 1) M80 
test aerosol, the low optical thickness r a (865) = 0.1, the low pigment concentration C = 0.1 
mg/m 3 and the geometry with 9 0 = 20°, 6 = 45.9° (Fig. 4.28); 2) U80 test aerosol, the 
high optical thickness r a (865) = 0.3, the high pigment concentration C — 1.0 mg/m 3 and 
the geometry 0o = 40°, 6 = 1° (Figure 4.29). We can observe that the fit of the computed 
total reflectance p*(A) — p r (A) to the given one was excellent, the objective function still 
being small (0.12%-0.40%) while the fits to the aerosol and water-leaving reflectances were 
slightly degraded. 

The retrieval error for the pigment concentration C was 6% for the positive calibration 
errors and about 9% for negative calibration errors. Apparently, the negative calibration 
errors generally affect the retrievals more adversely. The reason is the way we determine the 
Junge power-law parameter i/, i.e., using the Eq. (3.7), which results in the larger values for 
v when the errors are negative and the smaller when the errors are positive. The reflectivity 
of aerosols with the larger values of v increases and the total reflectivity is compensated by 
varying the water scattering parameter 6° in such a manner as to reduce the water-leaving 
reflectance keeping the pigment concentration almost unchanged. In contrast, the pigment 
concentration is varied when the v decreases. 

Compared to the results in the previous subsection, where the calibration errors were 
absent, the results were only slightly degraded. The average error of the retrieved pigment 
concentration for M80 case increased from 3.6% to 6% (the positive calibration error), for 
U80 case there was actually a decrease from 5.1% to 2.9% (the negative calibration error). 

From the above discussion we can conclude that the algorithm performance should be 
relatively insensitive to the radiometric calibration errors of the magnitude expected after 
in-orbit calibration adjustment. 



Normalized Water Leaving Reflectance, Normalized Water Leaving Reflectance, 

Aerosol and Combined Reflectance Aerosol and Combined Reflectance 



Figure 4.28 Retrieved the total, aerosol and water- 
leaving reflectances for the M80 test model with added 
positive calibration errors: r a (865) = 0.1, the pigment con- 
centration C = 0.1 mg/m 3 , the water scattering parameter 
6° = 0.3 m _1 , the viewing geometry 9q = 20°, 9 = 45.9° 



Figure 4.29 Retrieved the total, aerosol and water- 
leaving reflectances for the M80 test model with added 
positive calibration errors: r a (865) = 0.3, the pigment con- 
centration C = 1.0 mg/m 3 , the water scattering parameter 
6° = 0.3 m~ l , the viewing geometry 9q = 40°, 0 = 1° 
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RETRIEVAL OF OCEAN’S PIGMENT CONCENTRATION 

with positive calibration error 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M 80 , C 80 , T 80 , U 80 ). 

Also, provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 


, 

C; mg/m 3 

0.100 

0.500 

1.000 


result 


D 

result 

s d 

D 

result 

Sd 

D 

M80 

CDq = 0.993 

; 

x a (865) =0.10 

0.098 

1.6% 

2.00% 

0.494 



1.011 

4.2% 

1.10% 

x a (865) =0.20 

0.097 

2.2% 

3.00% 

0.487 

5.0% 

2.60% 

0.999 

8.5% 

0.10% 

x a (865) =0.30 

0.097 

2.2% 

3.00% 

0.485 

6.2% 

3.00% 

0.991 

11.% 

0.90% 



C; mg/m 3 

0.100 

0.500 

1.000 


result 

s d 

D 

result 

s d 

D 

result 

E9 

D 

C80 

C0 0 =0.988 

x a (865) =0.10 

0.098 

1.6% 

2.00% 

0.490 

3.6% 

2.00% 

0.995 

6.3% 

0.50% 

x a (865) =0.20 

0.097 

2.2% 

3.00% 

0.478 

7.1% 

4.40% 

0.966 

12.% 

3.40% 

x a (865) =0.30 

0.097 

2.6% 

3.00% 

0.471 

10.% 

5.80% 

0.940 

17.% 

6.00% 



C: mg/m 3 

0.100 

0.500 

1.000 


result 

Sd 

D 

result 

s d 

D 

result 

Sd 

D 

T80 

C0 0 = 0.953 

x a (865) =0.10 

0.095 

4.7% 

5.00% 

0.487 

2.4% 

2.60% 

1.037 

4.9% 

3.70% 

t a (865) =0.20 

0.093 

1.3% 

7.00% 

0.470 



0.935 

3.0% 

6.50% 

x a (865) =0.30 

0.096 

4.4% 

4.0% 

0.469 



0.943 

2.7% 

5.70% 



C: mg/m 3 

0.100 

0.500 

1.000 


result 

Sd 

D 

result 

E9 

D 

result 

Sd 

D 

U80 

COq =0.748 

x a (865) =0.10 

0.103 

0.8% 

3.00% 

0.519 

u% 

3.80% 

1.043 

2.3% 

4.30% 

x a (865) =0.20 

0.109 

2.8% 

9.00% 

0.543 

3.8% 

8.60% 

1.062 

2.9% 

6.20% 

x a (865) =0.30 

0.116 

5.7% 

16.0% 

0.573 

6.4% 

14.6% 

1.089 

3.6% 



Table 4.8 
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RETRIEVAL OF OCEAN’S PIGMENT CONCENTRATION 

with negative calibration error 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (M80, C80, T80, U80). 

Also, provided are the standard deviations over viewing geometries divided by the mean (S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/m 3 

0.100 

0.500 

1.000 


result 


D 

result 

a 

D 

result 

mm 

D 

M80 

COq =0.993 

x a (865) =0.10 

0.094 

2.3% 

6.00% 

0.479 



1.004 

0.8% 

0.40% 

x a (865) =0.20 

0.091 

4.0% 

9.00% 

0.463 

4.6% 


1.005 

2.2% 

0.50% 

r a (865) =0.30 

0.090 

4.7% 

3.00% 

0.452 



0.998 

4.0% 

0.20% 


. 

C: mg/m 3 

0.100 

0.500 

1.000 


result 


D 

result 


D 

result 

E9 

D 

C80 

C0 0 = 0.988 

x a (865) = 0.10 

0.094 

2.6% 

6.00% 

0.476 

2.5% 

4.80% 

0.997 

0.9% 

0.30% 

x a (865) =0.20 

0.090 

4.7% 

10.0% 

0.455 

5.2% 

9.00% 

0.983 

2.2% 

1.70% 

■c a (865) =0.30 

0.088 

5.3% 

12.0% 

0.438 

6.5% 

12.4% 

0.964 

4.4% 

3.60% 



C; mg/m 3 

0.100 

0.500 

1.000 


result 


D 

result 

mm 

D 

result 


D 

T80 

COq = 0.953 

x a (865) =0.10 

0.086 

10.% 

5.00% 

0.437 



1.038 

14.% 

3.80% 

x a (865) =0.20 

0.078 

14.% 

22.0% 

0.377 

19.% 

24.6% 

0.836 

3.4% 

16.4% 

x a (865) =0.30 

0.081 

12.% 

19.0% 

0.337 



KB9 

26.% 

36.0% 



C: mg/m 3 

0.100 

0.500 

1.000 


result 


D 

result 

19 

D 

result 


D 

U80 

C0 0 =0.748 

t a (865) = 0.10 

0.098 



0.496 

0.2% 

0.80% 

1.021 

2.7% 

2.10% 

x a (865) =0.20 

0.099 

2.7% 

1.00% 

0.507 



1.033 

3.5% 

3.30% 



3.2% 

4.0% 

0.530 

4.7% 

6.00% 

1.056 

4.1% 

5.60% 


Table 4.9 














































































































































D. Pseudo pA A) for U80 test aerosol models 
with vertical structure 
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So far it has been assumed that the aerosol is located within a thin layer just above the 
sea-surface. Actually, this is hardly ever the case in the natural conditions. In their study 
Ding and Gordon 43 came to a conclusion that as long as the aerosol was weakly absorbing 
zuq > 0.93, its vertical distribution had a negligible effect on the multiple scattering 
algorithm they used. On the other hand, if vj q decreased, the results degraded substantially. 
So, only highly absorbing aerosols need to be studied when the vertical atmospheric profile 
is taken into account. Later, Gordon 40 provided a simple explanation arguing that the 
increase of the amount of the Rayleigh gas in the aerosol layer increases the average path 
length of photons through the layer, and in this manner increases total absorption. In 
addition, the reflectance of the Rayleigh layer p r decreases for it is partially absorbed by 
the aerosol mixed up with the gas. 

While the aerosol vertical distribution is quite a complicated function, the easiest way 
to assess its impact on the radiative transfer processes is through the use of our two-layer 
atmospheric model where the top layer contains only the Rayleigh gas and the bottom layer 
is a mixture of the Rayleigh gas and aerosol with a constant mixing ratio. Denoting the 
fraction of the Rayleigh gas in the bottom layer as 7 , the optical thicknesses of the top and 
the bottom layers cam be simply expressed as: 

T.op = (1 ~ 7) T r, Tbot = 7 r r + T a , (4.3) 

where r r and r a are the total Rayleigh gas and aerosol optical thicknesses. The nomenclature 
for the U80 aerosol model is given in the Table 4.10 below. 

Table 4.10 Nomenclature for U80 models with vertical structure 


Notation 

Aerosol’s top altitude, km 

Fraction, 7 

U80 

up to 1 km 

0.0 

U280 

2.0 

UBM 

U480 

4.0 

mEm 
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The spectral variation of the aerosol reflectance p a + p ra for models U80, U280, U480 
with r a (865) = 0.2 is shown in the Figure 4.30. 



Wavelength, nm 

Figure 4.30 Spectral variation of p a +Pra for test models 
U80, U280, U480: r a (865) = 0.2, the viewing geometry 
0 O = 60°, 0 = 1° 


The aforementioned interpretation of this spectral behaviour can be clearly seen in this 
Figure. Indeed, the Rayleigh optical thickness r r in the blue is the largest comparatively to 
that in the NEEt, thus the aerosol has the lower reflectivity in the blue part of the spectrum. 
In the NIR, in contrast, the amount of Rayleigh scattering is very low and has a small 
impact on the reflectance field. 

To study the effect of the aerosol vertical distribution we generated data for the pseudo 
reflectance p t ( A) at the TOA using 7 values as indicated in the Table 4.10. Again, we 
applied the Solution 1. In building the system (3.10) only the first five bands were used. 
Particular results for the viewing geometry 9 0 = 20°, 9 = 1°, the optical thickness r a (865) 
and the ocean pigment concentration C = 0.5 mg/m 3 are shown in the Figures 4.11, 4.31, 




4.32 for two vertical distributions U80, U280, U480 respectively. 



Normalized Water Leaving Reflectance, Normalized Water Leaving Reflectance, 

Aerosol and Total Reflectance Aerosol and Total Reflectance 


0.040- 


U280, Q q = 20.0° 9 = l.Qff 0=90.0 P 


0.030-4 


Given Reflectances 

X 0 ° Retrieved Reflectances 



p,- P r 

c 

= 0.5 (mg/m 3 ) 

c’= 

0.604 (mg/rr? ) 




= 0.2 

*;= 

0.28 

O.Q25 - 



= 0.748 

©;« 

0.665 



b° 

= 0.30 

b °’= 

0.37 

0.020- 

Pra + Pa 




1.14% 


0.010-^ 





t Pw 


0.0 


Wavelength, nm 

Figure 4.31 Retrieved the total, aerosol and water- 
leaving reflectances for the U280 test model: r a (865) = 0.2, 
the pigment concentration C = 0.5 mg/m 3 , the water 
scattering parameter 6° = 0.3 m" 1 , the viewing geometry 
0 O =20°, 5 = 1° 



Figure 4.32 Retrieved the total, aerosol and water- 
leaving reflectances for the U480 test model: r a (865) = 0.2, 
the pigment concentration C = 0.5 mg/m 3 , the water 
scattering parameter 6° = 0.3 m“ l , the viewing geometry 
0 O = 20°, 9 = 1° 
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RETRIEVAL OF OCEAN’s PIGMENT CONCENTRATION 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and each of four hypothetical atmospheric aerosols (\I80, C80 , T80 , U80). 

Also , provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/ nr* 

0.200 

0.500 

2.000 


result 

KS 

D 

result 

s d 

D 

result 

s d 

D 

U80 

C0 0 = 0.748 

x a (865) =0.10 

0.100 

0 . 8 % 

0 . 00 % 

0.507 



1.028 

1 . 9 % 

2 . 80 % 

x a (865) =0.20 

0.104 

1 . 4 % 

4 . 00 % 

0.526 



1.053 

3 . 2 % 

5 . 30 % 

x a (865) =0.30 

0.110 

4 . 3 % 

10 . 0 % 

0.552 

5 . 8 % 

10 . 4 % 

1.070 

4 . 0 % 

7 . 00 % 


■V 

C: mg/m 3 

0.200 

0.500 

1.000 


result 

E9 

D 

result 

S d 

D 

result 

s d 

D 

U280 

C0 0 =0.748 

x a (865) =0.10 

0.110 

2 . 8 % 

10 . 0 % 

0.602 

9 . 8 % 


1305 

n .% 

30 . 5 % 

x a (865) =0.20 

0.122 

6 . 3 % 

22 . 0 % 

0.575 

3 . 0 % 

15 . 0 % 

1.120 

7 . 9 % 

12 . 0 % 

x a (865) =0.30 

0.138 

12 .% 

38 . 0 % 

0.631 

16 .% 

26 . 2 % 

0.959 

8 . 6 % 

4 . 10 % 



C: mg/m 3 

0.200 

0.500 

1.000 


result 

El 

D 

result 

El 

D 

result 

El 

D 

U480 

C0 0 = 0.748 

x a (865) =0.10 

0.120 

6 . 7 % 

20.0% 

0.941 

28 .% 

88.2% 

1.500 

0 . 0 % 

50 . 0 % 

t a (865) =0.20 

0.142 

14 .% 

42 . 0 % 

1.198 

19 .% 


1.500 

0 . 0 % 

50 . 0 % 

x a (865) =0.30 

0.172 

25 .% 

72 . 0 % 

1.084 

13 .% 

102.% 

1.500 

0 . 0 % 

50 . 0 % 


Table 4.11 


The detailed summary of the retrieval of the ocean pigment concentration is given in 
the Table 4.11 for three vertical distributions with optical thicknesses r a (865) = 0.1, 0.2, 0.3 
and three ocean pigment concentrations C = 0.1, 0.5, 1.0 mg/m 3 Note how the quality of 
retrievals degraded when the optical thickness r a increased. In fact, the standard deviation 
of the retrieval of the pigment concentration increased from 2.8% at r a (865) = 0.1 to 12% 
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at r a (865) = 0.3 for U280 test case and from 6.7% at r a (865) = 0.1 to 25% at r a (865) = 0.3 
for U480 test case (when C = 0.1 mg/m 3 ). 

The retrievals of the water-leaving reflectance had the deviation ranging from 10% 
to 30% for U280 and, from 20% up to 60% for U480 case. This is explained as follows. 
Increasing the layer thickness causes p a (A) + Pra(A) to become increasingly smaller as A 
increases. The algorithm deals with that through both choosing the largest value of m* 
available and increasing C to further reduce p t { A) - p r (A) in the blue relative to the red. 

The objective function S LS q on the average rose from 1.6% for U280 to 5.6% for 
U480 test case. No test case described in the previous subsections gave such large values 
for the objective function. This highlights the point that the vertical distribution for 
highly absorbing aerosols is a primary problem to deal with in any atmospheric correction 
algorithm. 

Encouraging is the fact that the retrievals for U280 test case, in general, produced 
acceptable results. This can be used to improve the algorithm’s performance when dealing 
with the aerosol vertical distribution. Indeed, if Look-Up Tables cure created with the Junge 
power-law models having a vertical distribution up to 2 km, it will have a little impact 
on the results when weakly absorbing aerosols are present. In contrast, when dealing with 
highly absorbing aerosols the results for cases when the aerosol is all near the surface and 
when the aerosol is mixed up with gas up to 4 km will be equally good. 
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E. Pseudo pA A) for Maritime-Mineral test aerosol models 

As has been noted before the Junge power-law aerosol models are characterized by a 
refractive index which does not vary with wavelength. To study the effect of a variable 
refractive index on our atmospheric correction results, we use the Mari time- Mineral test 
aerosol. Mineral component being of the crustal origin is frequently present in the maritime 
environment and is characterized by the trimodal log-normal size distribution. It is the min- 
eral component that is responsible for the variation of the refractive index with wavelength. 
Though, its fraction compared to the maritime component is rather small (from 5% to 20% 
being a strong function of relative humidity), its impact on aerosol radiative properties is 
significant. In Chapter 2 we have made a remark that the higher the modality of the size 
distribution, the better it is approximated by the Junge power-law size distribution. This 
fact was used in our simulations. 


The values of the refractive index for eight wavelengths and three Junge power-law 
parameters u are tabulated below. 


Table 4.12: Characteristics of the Maritime- Mineral test aerosol models. 




K9 

KrjjHQMEQQBI 

v = 2.0 

v = 3.0 

o 

II 

412 

1.53 — x 1.20e-2 

0.7679 

0.8990 

0.9297 

443 

1.53— t 9.10e— 3 

0.7978 

0.9185 

0.9445 

490 

1.53 — i 7.94e— 3 

0.8158 

0.9270 

0.9498 

510 

1.53 — i 7.34e— 3 

0.8252 

0.9315 

0.9528 

555 

1.53— i 5.40e— 3 

0.8576 

0.9469 

0.9637 

670 

1.53— i 4.30e— 3 

0.8863 

0.9566 

0.9687 

765 

1.53— i 3.16e— 3 

0.9146 

0.9670 

0.9752 

865 

1.53— i 1.20e— 3 

0.9645 

0.9866 

0.9898 


The spectral variation of r a (A)/r a (865) is almost independent of the refractive index 
and is similar to that shown in Fig. 4.4. The typical spectral variation of the reflectance 
p a (A) + /J ra (A) is plotted in the Figure 4.33 for the case with & 0 - 20°, 9 = \° ,<f> = 90° and 
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the optical thickness r a (865) = 0.2. We note immediately the behaviour of the Maritime- 
Mineral test aerosol with v = 2.0: 0 < e r < 1, cases that were not tabulated in the look-up 
tables (LUT), thus an extrapolation has to be performed in this case (remember that the 
selection of the set of {i/ y m r , m*} was based on the possibility of linear interpolation). Now, 
reminding that the calculation of u is based on the NIR bands exclusively (in the Solution 1), 
we might expect difficulties in solving the non-linear system of equations (3.4). Additionally, 
a strong dependence on the sun-viewing geometry might be expected (compare the results 
in Figures 4.34 - 4.36 to the results averaged over all sun-viewing geometries summarized 
in the Table 4.13) because of the increase in the ray path length. 



Figure 4.33 Reflectance p a ( A) + p ra (A) for the Maritime- 
Mineral test models with v = 2.0, 3.0, 4.0, r a (865) = 0.2, 
the viewing geometry 9 0 = 20°, 0 = 1° 


Encountered problems are reminiscent to those with the aerosol vertical structure dis- 
tribution from the previous subsection or the high reflectance T80 case from the Subsection 
B. Indeed, for the Maritime-Mineral case with v = 4.0 (see Fig. 4.34) the high reflectance 
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the algorithm achieves by increasing b° (which is constraint by the system (3.5)) and the 
aerosol optical thickness r a ( A) keeping C and too practically unchanged. This might not be 
the case, though, for other geometries. This leads sometimes to errors in C up to as much 
as 50% - 60%. 


RETRIEVAL OF OCEAN’S PIGMENT CONCENTRATION 

Tabulated are mean values of retrieved C for seven Sun-Viewing Geometries 
and hypothetical atmospheric Maritime-Mineral aerosol. 

Also, provided are the standard deviations over viewing geometries divided by the mean ( S d ) 
as well as the deviations from given parameters ( D ). 



C: mg/m 3 

0.100 

0.500 

1.000 
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E9 

D 
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D 

result 

WM 

D 
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mineral 

v = 2.0 

cOq = 0.965 

x a (865) =0.10 
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7.0% 

7.00% 

0.525 

5.8% 

5.00% 

1.018 

5.0% 

18.0% 

x a (865} =0.20 

0.115 

14.% 

15.0% 

0.577 

19.% 

15.4% 

1.074 

11.% 

7.40% 

x a (865) =0.30 
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22.% 

23.0% 
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1.256 
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25.6% 
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wm 

D 

result 

mm 

D 

result 

wm 

D 
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x a (865) =0.10 

0.110 

3.0% 

10.0% 

0.516 

2.2% 

3.20% 

0.959 

3.2% 

4.00% 

v = 3.0 

x a (865) = 0.20 

0.130 

11.% 

30.0% 

0.578 

9.2% 

15.6% 

0.948 

7.8% 

5.20% 

CD 0 =0.987 

x a (865) =0.30 

0.157 

24.% 

57.0% 

0.654 

11.% 

30.8% 

1.009 

17.% 

0.90% 



C; mg/m 3 

0.100 

0.500 

1.000 



result 

E9 

D 
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D 
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D 
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t a (865) =0.10 
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20.0% 

0.497 

1.4% 


0.914 

3.5% 

8.60% 

< 

ll 

o 

x a (865) = 0.20 

0.139 

14.% 

39.0% 

0.568 

11.% 

13.6% 

0.836 

4.8% 

16.4% 

CD 0 =0.990 

x a (865) =0.30 

0.195 

34.% 

95.0% 
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Table 4.13 



















































































































Normalized Water Leaving Reflectance, Normalized Water Leaving Reflectance, 

Aerosol and Combined Reflectance Aerosol and Combined Reflectance 



Wavelength, nm 

Figure 4.34 Retrieved the total, aerosol and water- 
leaving reflect ances for the Maritime-Mineral test model 
with v = 2.0: r 0 (865) = 0.2, the pigment concentration 
C — 0.5 mg/m 3 , the water scattering parameter 6° = 0.3 
m -1 , the viewing geometry 9 0 = 20°, 0 = 1° 



Figure 4.35 Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime-Mineral test model 
with v = 3.0: r a (865) = 0.2, the pigment concentration 
C = 0.5 mg/m 3 , the water scattering parameter 6° = 0.3 
m~ l , the viewing geometry 0q = 20°, 0 = 1° 
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Wavelength, nm 

Figure 4.36 Retrieved the total, aerosol and water- 
leaving reflectances for the Maritime- Mineral test model 
with v — 4.0: r 0 (865) = 0.2, the pigment concentration 
C = 0.5 mg/m 3 , the water scattering parameter 6° = 0.3 
m -1 , the viewing geometry &o = 20°, 0 = 1° 


While the values of the single-scattering albedo rcr 0 (865) were within 5% - 6% range in 
the NIR, they made no sense in the blue, particularly for v = 2.0 test model. 

The objective function was on the average equal to 1.2% (again, in building the system 
(3.10) only the information from the first five bands was used). In some individual cases 
•Slsq was as small as 0.4 — 0.8% making it difficult to use this objective function for an 
identification of the maritime-mineral aerosols. Note, that with the Junge power-law models 
we cannot estimate the spectral variation of the refractive index, so just observing a large 
objective function value 5 lsq will not recognize whether it happened because of the aerosol 
vertical distribution or because of the spectral variation of the refractive index. 




V. TESTING SPECTRUM OPTIMIZATION WITH SeaWiFS 


The SeaWiFS, the Sea-viewing Wide Field-of-view Sensor, was launched by NASA in 
August 1997. Its purpose is to examine oceanic factors that affect global change and to 
assess the ocean’s role in the global carbon cycle, as well as other biogeochemical cycles. 8 

Major instrument parameters have been summarized in the Table 5.1. The selection of 
of the SeaWiFS bands was based on the spectral absorption characteristics of some co mm on 
in- water optical constituents, mean extraterrestrial irradiances, and spectral transmittance 
of the atmospheric constituents oxygen and water vapor. SeaWiFS produces two types of 
science data: LAC (Local Area Coverage), with a nominal ground resolution of 1.13 km 
(1,285 picture elements, or pixels, per scan line) and GAC (Global Area Coverage), with a 
nominal ground resolution of 4.5 km, which is subsampled from LAC data every fourth line 
and every fourth pixel (such that is comprised of 248 pixels per scan line). There are four 
Levels of SeaWiFS data products: 

0 Raw radiance counts from all bands amd instrument telemetry. 

1A All Level-0 data, appended calibration and navigation data, and reformatted instru- 
ment and selected spacecraft telemetry. 

2 Derived from corresponding Level-IA GAC product by applying the sensor calibration, 
the atmospheric corrections, and bio-optical algorithms. 

3 Binned (a global, equal-area grid whose grid cells, or “bins”, are approximately 81 
km 2 ) data product that consists of the accumulated data for all Level-2 GAC products 
corresponding to a period of one day, 8 days, a calendar month, or a calendar year. 
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Table 5.1: Major Parameters and Characteristics of the SeaWiFS ocean color sensor. 


Instrument Bands 

Band 

Center 

Wavelength 

Saturation 

Input 

SNR 


Wavelength 

FWHM 

Radiance 

Radiance 



[nm] 

[nm] 

[mW cm~ 2 /zm~ 1 sr _1 ] 

[mW cm'^m' 1 sr _1 ] 


1 

412 

402-422 

13.63 

9.10 

499 

2 

443 

433-453 

13.25 

8.41 

674 

3 

490 

480-500 

10.50 

6.56 

667 

4 

510 

500-520 

9.08 

5.64 

640 

5 

555 

545-565 

7.44 

4.57 

596 

6 

670 

660-680 

4.20 

2.46 

442 

7 

765 

745-785 

3.00 

1.61 

455 

8 

865 

845-885 

2.13 

1.09 

467 

Sensor Accuracy 

Radiance Accuracy: 


<5% absolute each band 


Relative Precision: 


<1% linearity 



Between Band Precision: 

<5% relative band-to-band (over 0.5- 0.9 full scale) 

Polarization: 


<2% sensitivity (all angles) 


Nadir Resolution: 


1.1km LAC; 4.5km GAC 


Mission ChcLracteristics 

Orbit Type: 


Sun Synchronous at 705 km 


Equator Crossing: 


Noon ±20 min., descending 


Saturation Recovery: 

<10 samples 



Modulation Transfer Function: 

>0.3 at Nyquist 



Swath width (at equator): 

2,801 km LAC (±58.3°) 





1,502 km GAC (±45.0 <> ) 


Scan Plane Tilt: 


4-20°, 0°, —20° 



Digitization: 


10 bits 




It is the Level-2 data product generation, or more exactly, its atmospheric correction part, 
that we will be concerned with in this chapter. 

To facilitate a rapid processing, display, analysis, and quality control of all SeaWiFS 
products, GSFC NASA created a comprehensive image analysis product named SeaDAS . 
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The NASA Standard Atmospheric Correction algorithm has been developed by Gordon 
and Wang. 11 It will not be described here since it can be found in published literature 
elsewhere. 11 ’ 45 This algorithm has proven to be fast, reliable and simple to implement. 
Our purpose will be to compare the performances of the aforementioned NASA Standard 
Correction and the Spectrum Optimization Algorithm (described in Chapter 2) within the 
SeaDAS image processing system. 


Table 5.2: 

Comparative Methodology between the 
NASA Standard and Spectrum Optimization Algorithms. 



Method or 

Structure 

Function 

NASA Standard Correction 

Spectrum Optimization 

Aerosol 

Assumptions* 

12 models 

for M,C,T models at 4 relative 
humidities 

72 models 

for Junge Power-Law models: 
6 values of u, 2 of m r , 6 of m* 

Data Storage! 

LUT for p*(p*.) + pra(p*,) 

LUT for Pa(r a ) + Pra(r a ) 

Aerosol 

Properties: 

Bands 7 and 8 are used to determine 
e* (765, 865) 

LUTs are used to bracket the best 
2 models. Further interpolation pro- 
vides an “average” model. 

Bands 7 and 8 are used to determine 
v : 

given ?'(765, 865). Then r a (865) 
is determined for each m r ,mi combina- 
tion using Newton-Raphson iteration. 

Aerosol 

Reflectance 

Pra(Aj) 

Based on extrapolation of the 
“average” model into the visible 

Quadratic interpolation for 
scattering parameters for LUT models. 

Water-Leaving 
Reflect. Pu»(A) 

Computed directly 

Non-linear optimization using 
Bands 1 through 5 (Band 6 is optional) 

• 

Diffuse 

Transmit. t(A) 

l 

Bidirectional effects not included 

Bidirectional effects not included 


* M, C, T, U shortcuts stand for Maritime, Coastal, Tropospheric, Urban models 
t p<x» stands for the single-scattering reflectance. 
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From the Table 5.2 we can make a note on the way the pigment concentrations are 
obtained by each of the algorithms: the NASA Standard Correction algorithms utilizes two 
bands (centered at 490 nm and 555 nm), while the Spectrum Optimization algorithm uses 
the first five bands (the usage of Band 6 is optional). This is why we might expect better 
retrievals of the pigment concentration with the latter algorithm. 

For atmospheric correction experiments we select two L1A-GAC images that were taken 
on October 6, 1997 (corresponds to the 279th day of a year) and October 8, 1997 (DAY-281). 
We consider the same grid cell with its top left corner at (Lat,Lon) = (42°84”,78°43”) and 
the bottom right comer at (29°12 n ,63°77”) Since the time difference between the pictures 
is only two days and, the “productivity” (i.e. the production of organic material by algae 
and some bacteria that exist at the lowest levels of the oceanic food chain that utilize the 
sunlight or chemical energy as their source of energy) during fall is not high, the change 
in the pigment concentration is not expected to be large. DAY-279 is characterized by a 
highly turbid atmosphere, in contrast to the clear sky DAY-281. The true color pictures of 
those days are shown in Fig. 5.1. Figures 5.2 and 5.3 show Level-IA raw data for each of 
eight bands. After applying calibration and other metereological data, we perform Level-2 
product generation using both algorithms. 

The results of retrieving the pigment concentration are demonstrated in Figures 5.4 - 
5.7. It can be noticed that the pigment concentration for DAY-281 obtained by NASA 
Standard Correction and Spectrum Optimization algorithms differ very little. Indeed, 
during DAY-281 we expect the aerosol to be completely of Maritime nature and of very 
small optical thickness, the case which both algorithms performed very well during initial 
testing. A sharp contrast in performance appears while processing the DAY-279 image. 
From Figure 5.5 we note a drastic change in pigment concentration comparatively to that 
of DAY-281 produced by NASA Standard Correction Algorithm. More precisely, retrieved 
concentrations for DAY-279 were from 2 to 3 times as large as that for DAY-281. Spectrum 
Optimization Algorithm results (Figures 5.6, 5.7) are substantially better, showing no 
large variations in the pigment concentration for DAY-279 and DAY-281. Needless to 
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say, Spectrum Optimization results are typically much smoother comparatively to results 
obtained by NASA Standard Correction Algorithm. 

We further analyzed the water-leaving radiances retrieved by the algorithms. In or- 
der to escape the contamination of pictures by clouds, we extract subimages from our 
L2-GAC products. Histogram analysis is performed on water-leaving radiances for Band 
2 (443 nm) and 5 (555 nm) for both days (the top corner of the subimage is located at 
(Lat,Lon)=(39°60” ,70°33”) and the bottom right corner at (Lat,Lon)=(37°80”,67°32”). 

NASA Standard Correction Algorithm: 

For DAY-281 this algorithm (Fig. 5.8b) indicated two peaks of the water-leaving 
radiance for Band 2: first at 1100 (mW cm -2 um -1 sr” 1 ) with the full width at half 
maximum (FWHM) of about 200 (mW cm" 2 um" 1 sr" 1 ) and the second at 1390 (mW cm" 2 
um" 1 sr" -1 ) with FWHM of about 110 (mW cm" 2 run -1 sr” 1 ). For DAY-279 FWHM’s were 
much wider (the second peak was almost unidentifiable) and shifted about 500 (mW cm"" 2 
um" 1 sr"" 1 ) toward higher values of radiances, suggesting that the atmospheric correction 
algorithm underestimated the aerosol component. 

For Band 5, both days showed a single peak, although that for DAY-279 was again 
shifted toward higher values and had FWHM more then triple that of DAY-281. 

Spectrum Optimization Algorithm: 

For DAY-281 this algorithm (Fig. 5.10b) for both bands indicated the results similar to 
those obtained by NASA Standard Correction Algorithm. The results differ for DAY-279. 
FWHM were almost unchanged, though the radiances shifted toward higher values as well. 
The second radiance peak for Band 2 is clearly visible (the noise in the histogram may be 
the result of the image’ small size). 
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SI 997279 17 1919 . L1A_HNSG_BRS S19972 8117 03 57 . L1A_HNSG_BRS 

DAY - 279, October 06, 1997 DAY - 281, October 08, 1997 

Figure 5.1 “True color” images taken by SeaWiFS on October 6, 1997 (Day-279) 
and October 8, 1997 (Day-281). 
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Figure 5.2 Raw data images for Bands 1 through 8 taken by SeaWiFS on 
October 6, 1997 (Day-279). 












Figure 5.3 Raw data images for Bands 1 through 8 taken by SeaWiFS on 
October 8, 1997 (Day-281). 
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Figure 5.4 Level-2 pigment concentration obtained using the NASA Standard Al- 
gorithm for Day-281 
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Figure 5.5 Level-2 pigment concentration obtained using the NASA Standard Al- 
gorithm for Day-279 
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Figure 5.6 Level-2 pigment concentration obtained using the Spectrum Optimization 
Algorithm for Day-281 


Figure 5.7 Level-2 pigment concentration obtained using the Spectrum Optimization 
Algorithm for Day-279 
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Figure 5.8 Level-2 Band 2 (443 nm) image obtained using the NASA 
Standard Correction Algorithm for a) Day-279, b) Day-281 and, corre- 
sponding histograms 
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Figure 5.9 Level-2 Band 5 (555 nm) image obtained using the NASA 
Standard Correction Algorithm for a) Day-279, b) Day-281 and, corre- 
sponding histograms 
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Figure 5.10 Level-2 Band 2 (443 nm) image obtained using the Spec- 
trum Optimization Algorithm for a) Day-279, b) Day-281 and, corre- 
sponding histograms 
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Figure 5.11 Level-2 Band 5 (555 run) image obtained using the Spec- 
trum Optimization Algorithm for a) Day-279, b) Day-281 and, corre- 
sponding histograms 
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Conclusion . The results presented in this chapter show that: 

- Spectrum Optimization algorithm (SO) was successfully implemented in an image 
processing environment. 

- SO produced excellent C for both days (in spite of turbid atmosphere on DAY-279) in 
agreement with the predictions in Chapter 4. 

- SO produced normalized water-leaving radiances (nLw) with narrower distributions 
than NASA Standard algorithm and with smaller shift to higher nLw on DAY-279. 

- compared to the results of NASA Standard Correction Algorithm, SO produced less 
noisy images, i.e., more water features were observed (the aspect is not understood at 
present). 



VI. CONCLUSION 


An atmospheric correction algorithm based on the non-linear optimization procedures 
capable of retrieving simultaneously the aerosol and ocean water optical properties under 
the presence of both highly and weakly absorbing aerosols has been developed in this 
dissertation. The Junge power-law with variable refractive index have been choosen as the 
model aerosols, since their radiative properties can be varied continuously via the parameters 
{z/, m — providing a possibility to make use of non-linear optimization techniques. 

The algorithm performance with the simulated data revealed the robustness of the 
method when dealing with both highly and weakly absorbing atmospheric aerosols. In fact, 
the objective function 5 LS q, the “measure of the quality of the data fit”, was extremely 
small, i.e., 0.27% for the Maritime-80 and 0.15% for the Urban-80 test aerosol models. The 
retrieval of the ocean pigment concentration was surprisingly accurate, the error reaching 
on the average 3.5% for Maritime-80, and 5.1% for Urban-80 models. The error in the 
single-scattering albedo zuq was on the average 3%. As long as the phase functions of the 
Junge power-law aerosol models does not reproduce properly the phase functions of realistic 
aerosols, the retrieval of the optical thickness r(A) was noticably poorer then that of other 
parameters. The error was sometimes as high as 30-40%. Unfortunately, this is a drawback 
of the algorit hm . On the other hand, the fact that the other parameters, such as the 
pigment concentration and the single-scattering albedo, were obtained within an acceptable 
accuracy sets forward a suggestion as to whether the actual form of the phase function is 
necessary for that purpose. 

We studied further the effects of the calibration errors after effecting an in-orbit cali- 
bration adjustment. The result of retrievals was not substantially worse. In the case of the 
pigment concentration the error was on the average 6% for positive calibration errors and 
9% for the negative ones. 
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The presence of the aerosol vertical structure has adversely affected the results of the 
atmospheric corrections. The impact was so strong in cases when the aerosol was uniformly 
mixed with the Rayleigh gas that the results could have hardly been called valid. It is 
suggested, however, to use the Junge power-law models which are partially mixed up with 
gas. In this case, the error of retrievals is expected to be that of the error when calibration 
uncertainty is present in the radiance data. 

The effects of the variable refractive index with respect to the wavelength has been 
also analyzed. The maritime-mineral aerosols are appropriate representatives of this class 
of aerosols. The study shows that the retrieval of the pigment concentration is still accept- 
able (the error was within 25%), though the accuracy of the other parameters degraded 
substantially. However, both cases, the aerosol vertical distribution and the refractive in- 
dex variable with the wavelength were easily recognized by the algorithm. The objective 
function 5 LS q was on the average 10-15% for the aerosol vertical distribution and 1.2% for 
the maritime-mineral aerosols. 

The algorithm has been successfully incorporated into the SeaWiFS image processing 
system SeaDAS. Images of the sea surface have been corrected for two days with two 
distinctive aerosol present: turbid and non-turbid. In both cases, the retrieval of the 
pigment concentration produced similar results, in contrast, to those obtained with the 
NASA Standard Correction Algorithm, which failed in the case of the turbid atmosphere. 
The water-leaving radiances for two days were compared with the help of image histograms 
that revealed a general shift of radiances’ peaks in all bands toward higher values which 
might be a consequence of an incorrect calibration adjustment supplied to the algorithm. 

Though the extensions to the algorithm are yet to be studied, this work demonstrates 
a general applicability of the method to ocean color remote sensing. 



APPENDIX: Foundations of Quasi-Newton Optimization Metkods 


In this Appendix laid out will be the basic theoretical foundation of the non-linear pro- 
gramming (optimization) methods belonging to the so-called quasi-Newton class of methods. 
The theory will be supported with an examplary Davidon-Fletcher-Powell algorithm solu- 
tion. Apparently, it is impossible to describe in a small section all the cornerstones of this 
wide field of computer science. So, the purpose is limited to highlighting the major building 
blocks of any optimization procedure, the knowledge of which provides the pavement for 
better and ellegant problem understanding and solving. 


A.l Notation, Terminology and Definitions. 


Consider the n- dimensional euclidian space E n . A point x = [xi , . . . , x n ] T in this 
space is a column vector of its components xi, . . .,x n . In general, it will be assumed 
that the problem arising functions f (x) are smooth, that is continuous and continuously 
differentiable (C 1 ). Let us have two other C 1 functions hj(x) and ^(x). Then, the non-linear 
optimization problem might be stated as follows. 

Minimize 

/ (x), x £ E n (A.l) 

subject to m linear and/or non-linear constraints: 

/ij(x) = 0, j = 1, ... ,m (A2) 

and ( p — m) linear and/or non-linear non-equality constraints: 

< 7 j(x) > 0, j = m + l,...,p (-4-3) 
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Now, a function <f>(x) is called convex over the domain 7Z if for any two vectors, say xi 
and x 2 6 1Z the following inequality is satisfied: 

<f>(dx x + (1 - 0)x 2 ) < ^( Xl ) + (1 - 0)0(x 2 ), 0 < 9 < 1 (A.4) 

In other words, a convex function 0(x), by its classical definition (A.4), cannot have any 
value larger then the values of the function obtained by linear interpolation between <£(xi) 
and <^(x 2 ). A differentiable convex function has the following properties: 

a ) ~ ^(*i) > V T ^(xi) (x 2 - xi) for Vxi,xi 

b) the Hessian matrix (see below) is positive definite 

c) ^(x) has only one extremum. 

The convexity of a set of points {x^} is defined the same way, i.e. all the vectors 
x = 9x x + (1 — 0)x 2 , Vxi,x 2 , 0 < 9 < 1, should belong to the set {xj}. Other seven types 
of convexity, very relevant to the optimization problems, have been discussed by Ponstein. 46 

For the norm of a vector ||x|| we choose the common i 2 - norm, defined as 

Mi = 

Note, as has been pointed out by Powell, 47 ii-norm, defined as ||x||i = x i* vei T often 

leads to bad- conditioning in the approximation methods. 

Iterative procedures generate a sequence of points x^^x^ 2 ), . . . ,x ^ k \ . . ., or {x^} 
(where k indicates the iteration number) converging to a fixed, 48,49 stationary point x* 
which provides the solution (or optimal value solution) to a problem. 35 

For any iteration C l -function /(x^) at the iteration point x^ there is a vector of its 
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first partial derivatives, or gradient vector, V/(x^ fc ^) 


df(xW)/d Xl 

df(x^)/dx 2 

df{xW)/dx n ) 


= V/(x( fc )) = V/ (fc > 


(A.6) 


If /( x is twice continuously differentiable (C 2 ), then there exists a matrix of its 
second partial derivatives or, Hessian matrix H(x^), written as 


df(xW)/dz{ ... df(xW)/d Xl dx n ' 
[ df(x ( - k) )/dx n dxi ... df{x^)/dx 2 n , 


= V 2 /(x (fe >) = H(x (fc >) (A. 7) 


This matrix is square and symmetric. In the n- dimensional Euclidian space we define a line 
as a set of points 


x = x(A) = x' + As 


(A.8) 


for all A, a fixed point x! and any vector s collinear with, the line (see Fig.A.l). 


From (A.8) we can determine derivatives of /(x) along the line: 


d r-^ dxA A) d d t_. 

Tx = X’*T97 { =X'“d^ = s V 


(A.9) 


The slope of /(x) along the line at any point x(A) is 


d A 


= s T V/ = V T / s 


(A.10) 


The curvature along the line is 




(A. 11) 
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where V/ and V 2 / are evaluated at x(A). 



Figure A.l Definition of a line. 

Assume that is some point of an iteration. Expanding /(x) in the vicinity of x( fc ) 
into the Taylor series, truncating all terms of orders higher then two, we obtain: 

/(x) - /(x( fe >) a (V/(x^)) T (x - x (fe) ) + i (x- x <*> f v 2 /(x< fc >) (x - x< fc >) (A. 12) 

where V 2 /(x^ fe ^) = H(x^ fc ^) is the Hessian matrix. Note, for a quadratic unconstrained 
functions (Eq. (A. 12)), the necessary conditions for x* to be a local minimum is 

1. /(x) is differentiable at x*, 

2. V/(x*) = 0, 

3. ( sufficient condition ) the Hessian is positive definite. 


A.2 Steepest Descent Method. 

Before continuing, it can be observed that the iteration step from the point x^ fc ) to 
may be viewed as 

x (fc+i) _ X W + Ax (*) _ x (*> + \W s( fe ) = + \’( k) s (fc) (A. 13) 

where s^ k ^ is a uni t vector in the direction of Ax^ k \ s^ k ^ is a generic vector in the direction 
Ax( fc \ \( k ^ and \ m ( k ) are scalars such that Ax^ fc ^ = A^ = A*( fe ) s( fe ). 


In the search for a reliable direction one may suggest to pick it up in the direction 
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of the greatest decrease of /(x) with respect to its variable x, i.e. 


z(») ■ V /(* (fe) ) 

nv/(xW)ir 


(A. 14) 


where ||V/(x (fe ))|| is the chosen norm of a vector. Plugging this back into Eq. (A. 13), we 
obtain the following step 


x (fc+i) _ x (fc) _ A (fc) 


V/(x^) 

l|V/(x(*))|| 


= x (fe) - \< k) V/(x< fe >) 


(A.15) 


where A*( fe ) is a variable parameter in the Eq. (A. 8). In order to provide the best step 
possible, it is needed to perform a one- dimensional (in A*( fc )) search for a minimum of 
/( X ( fc+1 >) (the concept of the one-dimensional search has been well-developed and can 
be found elsewhere in literature. However, in the author’s opinion, the one described by 
Armijo 50 proved to be very fast and extremely robust in real problems). Thus, we arrive 
at the point x^ fc+1 ^. Proceeding this way with x^ +2 ^, . . . , the optimal value solution x* 
may be found that satisfies some requirements for a termination: value of /(x), x itself, 
smallness of A, zeroing of V/(x), etc. At the solution point x*, typically, the Hessian is 
determined and, if it is not positive definite, the x should be moved from its position one 
way or another (or penalized), like in the case of having a saddle point. 


Returning back to the Eq. (A.15), we may observe that formally A can be found 
from the condition of minimum of a function (note, this is not generally true for constrained 
problems): 


d/(x( fe > + A( fe ) s< fe >) 
d A 


= 0 


(A.16) 


If /(x) is a quadratic function (i.e. representable by (A. 12)), then from (A. 12) 


(V/(xW)) T s< fe > + (sW) r H(xW) A s<*> = 0 (A.17) 


So that, 

(V/(x< fe >)) r s< fe > 

“ (s( fc )) T H(x(*>) s( fe ) 


(A.18) 
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A. 3 Newton Method 

Let us try to look at the minimum of /(x) in the direction of i.e. 


d /(x< fe >) 
d Ax 

(A.19) 

From (A. 12) this yields 


0 = V/(x (fc >) + V 2 /(x (fc) ) x (fc) 

(A. 20) 

so that, 

J» _ V/(x<*>) 

V 2 /(xW) 

(A. 21) 

or, 


x (fe+1) = x (fe) - [V 2 /(x (fc) )] -1 V/(x (fc) ) 

(A.22) 


This would provide the solution if only /(x) were strictly quadratic. Generally speaking, it 
is not ever the case. Conforming (A. 22) to (A.15), we obtain the following 


M) = m _ xW [V7(xW)]-* [V/(x( fe >)] 

|| [V*/(arf*))]- 1 [V/(x(^))]|| 


(A.23) 


or, 

x (fe+1 > = x (fe) - H- x (x (A: >) V/(x< fc >) (A.24) 


Observe that in contrast to the Steepest Descent method the search direction is now given 

by 

s (fe) = -H _1 (x (fe) ) V/(x (fe) ) (A.25) 


Well, numerical complexities of this approach are apparent: 1) the second derivatives should 
be computed; 2) a matrix inversion is required. All this makes the method extremely 
impractical. 
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A. 4 Quasi-Newton or Variable Metric Methods 

This section will provide the way to approximate the inverse Hessian matrix H -1 (x^ fe ^) 
without actually computing it (from here the term “quasi-Newton”). 35,38 Denoting 

«(x<*>) « (A.26) 

the Eq. (A. 24) becomes 

x (fc+1) = x (fe) - A* (fc) $(x (fe) ) V/(x (fc) ) (A.27) 

(note, $(x^ fe ^) corresponds to the identity matrix in the Steepest Descent method and, is 
typically called the direction matrix). 

From analytical geometry it is known that for quadratic functions the Hessian matrix 
is sy mm etric. Using (A.12) and noting that H(x^ fe ^) is set, i.e. constant, and, constant also 
are the (V/(x^ fc ^)) T and /(x^ fe ^), we find the gradient of /(x) to be 

V/(xf= (V/(x< fe >)) T + \ (H(x< fc >) + H r (x< fe >)) (x - x<*>) 

mM 

= (V/(x< fe >)) r + H(x (fc) ) (x - x (fe) ) (A.28) 

where we used the rule of the gradient of the scalar product of vectors x and Hx. 

Premultiplying both sides by H -1 (x^ fe ^), 

x (fc+1) - x< fe > = H- X (x (fe )) [V/(x< fe+1 >) - V/(x< fe) )] (A.29) 

Trying to find A$( fe ) 

H- X (x (fe) ) w u> $ (fe+1) = w ($ (fc) + A$ (fc) ) (A. 30) 

where a; is a scaling parameter, it is introduced since functions are typically non- quadratic, 
the technique we already used before in (A. 24). At stage ( k + 1) we are provided with 
x< fe \ V/(x< fe >), V/(x( fe+1) ) and Let Ag< fe > = V/(x< fe+1) ) - V/(x< fc >). Then 

Ag^ fe ^ = — Ax^ 


(A- 31) 
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or, 


A $ (<!) Ag (fe) = iAx (fc) - $(*) Ae (fe) 

u> s 


(A. 32) 


Eq. (A. 32) may be viewed as a system of n 2 equations in n 2 variables. Once it is solved, 
the information on A$^ fe ) becomes available. By direct substitution it can be verified that 


= I Ax(fe) y T _ Ag(*> 2 r 

w y T Ag( fe ) z T Ag( fe ) 


(A. 33) 


is the solution where y and z are arbitrary vectors. The freedom of choosing y and z leads 
to a variety of quasi-Newton methods. Some examples are: 


Brovden : 


y = z = Ax^ — Ag^ 


A$ 


(*) 

Broyden 


[Ax^) — Ag( fc )] (Ax^))^ 1 
Ag( fe ) (*(*)) T xW 


(A. 34) 
(A. 35) 


Davidon-Fletcher-Powell' 


A¥ (fc) 


y = Ax (fc) , z = $( fc )Ag (fe ) 

[Ax^] [Ax( fe )] r [Ag<*>] [Ag^)] 2, [$( fc )] T 

DFP “ [Ax( fe )] r [Ag( fe >] [Ag( fc )] T $( fe > [Ag<*>] 


(A.36) 

(A.37) 


Brovden- Fletcher- Goldfarb- Shanno : 

( [AgW] T $( fe ) Ag( fc > Ax<*> [Ax< fc >]‘ ] 
l 1+ Ax( fc >Ag( fe ) [Ax( fc )] r AgW J 

f x (*) [Ag( fe )] T + $( fc ) [Ag( fe ^] [x( fe )] r 
{ [x< fe )p Ag<*> 


(A.38) 


A$( fc ) 


BFGS 




It should be noted, however, that BFGS algorithm has been derived from other suggestions 
and is given here for completeness since it is complement ary to DFP algorithm. 

Once the A$( fc ) are computed, (A. 27) should be used to determine the length of the 
”next” step Ax^ through one- dimensional search, namely, minimizing \*( k \ Again, the 
iteration procedure is repeated until some of terminating conditions are met. 
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A. 5 Example: analytical DFP minimization 

As an example let us try to minimize an unconstrained function of two variables: 

/(* i , x 2 ) = 2 (ii - l) 2 + 3 ( x 2 - 2) 2 (Ae.l) 

O 

This simple quadratic function reaches its minimum f(xi, x 2 ) = 0 at the point x* = 
[1 2] t . We are trying to use the information from the previous subsections to perform 

all the required iterations assuming that the function and its derivatives can be computed 
analytically everywhere. 



and, 


- - 0 = -580 + 2968 A* (0) 

d A*(°) 


(Ae.6) 
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So that, A*(°) = 0.195418. From (Ae.3) we determine the first iteration point x^ 1 ); 


x* 1 ) = * 

- ,(D 

X 2 


1.873315 

1.482476 


We complete the step with the other relevant results: 


/(x< l >) = 2.3288 Ax<°> = [ 3 126685 

V ' [—3.517524] 

Ag(°) = [~ 12 - 506740 ] $ (o) 

8 [-21.105144] 


-3.1266851 (1) [ 3.493260" 

-3.517524 n ’ -3.105144 


The update to the approximate inverse hessian matrix can be found using (A.37): 

[-3.126685 0] [-3.126685 -3.517524] 

A *(0) = [-3-517524 Oj [ 0 0 J 


-3.126685 - 3.517524] 


1 0] [-12.506740 0] [-12.506740 

0 1] [-21.105144 0] [ 0 

[-12.506740 -21.105144] J J 


-21.105144 

0 

~— 12.506740 
-21.105144 


IM 

40' 

44 


-0.173645 -0.341543 
-0.341543 -0.630937 


And, 


[ 0 . 
[- 0 . 


.826355 -0.341543 

.341543 -0.369063 


Step 2. 


As before, the formula for the step function (A.27): 


x< 2 > = 


x[ 2) ' [1.873315] _ (1) 

x[ 2) [l .482476] 


3.947213" 

-2.339044 


The step size: 


/(x (2) ) = 2.328851 - 21.051708 A* (1) + 47.574360 (A* (1) ) 2 


(Ae.ll) 
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and, 


d /(x f2) ) 

d \< l ) 


= 0 = -21.051708 + 95.148720 A*< 1 > 


(Ae.12) 


x 2 



-4 -2 0 2 4 6 


Figure A. 2 Steps in minimization of function f(x i , X 2 ) = 
2 (ri — l) 2 + 3 ( 2:2 — 2) 2 using DFP-algorithm. 
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Thus, 


x (3) = [1.00 2.00] /(x (3) ) = 6.71 x 1CT 15 


Discussion . 

Even this simple example shows some general behaviour of DFP (and other quasi-Newton) 
iterations. The function values become minimized much faster well before the iteration point 
reaches the desired optimal value point. The fact that in this demonstration we reached the 
minimum in three steps is the result of the given function /(x) being quadratic. Indeed, 
quasi-Newton methods on each iteration select the direction of the appropriate conjugate 
directions (similar to Conjugate Gradient methods that were not considered here) and, in 
case of two-dimensional function there are just two of such directions. 
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Abstract 


A radiometric system, deployed from a ship, is used to measure directly the influence of the 
presence of breaking waves (whitecaps) on the upwelling radiance above the sea surface. 
Estimates of their remote-sensing augmented spectral reflectance, i.e., the temporally-averaged or 
spatially-averaged increase in the ocean's reflectance over and above the reflectance in the 
absence of breaking waves, are provided from measurements in the Tropical Pacific. The 
accuracy of these estimates is dependent on the ability to determine radiometrically the 
background reflectance of the water. In the visible, the remote-sensing augmented spectral 
reflectance of whitecaps measured in the open ocean is found to be essentially independent of 
wavelength, -0.001 at a wind speed of -9 m/s, and -0.002 at -12 m/s. This is a factor of -2.5 
less than predicted earlier based on the statistical relationship between fractional coverage and 
wind speed, and the estimated average reflectance of individual whitecaps. In the near infrared 
(860 nm), the remote-sensing augmented spectral reflectance falls to 80% of its value in the 
visible. 


Introduction 

Due to the relatively small radiance backscattered into the atmosphere from below the ocean 
surface (the water-leaving radiance) in comparison to the total light arriving at a space-borne 
ocean-color sensor, the usefulness of ocean color data depends primarily on the accuracy of 
atmospheric correction algorithms [ Gordon and Morel, 1983]. Atmospheric correction over the 
ocean is effected by measurement of the upwelling signal arriving at the satellite at spectral bands 
where the radiance backscattered out of the ocean is known to contribute very little or no signal 
[Gordon and Wang, 1994a]. Utilization of spectral bands beyond -740 nm in Case 1 waters, i.e., 
waters in which the optical properties are determined by phytoplankton, their decay products, and 
the water itself [Gordon and Morel, 1983], can be justified in the absence of whitecaps due to the 
strong absorption properties of liquid water in the near infrared (NIR). When waves break on the 
sea surface, the presence of the resulting whitecaps occurring within a satellite image pixel 
increases the radiometric signal by an amount dependent on the whitecap coverage and their 
reflectance. This increase in apparent surface reflectance over and above that in the absence of 
breaking waves is referred to here as the augmented reflectance. [In this paper, we consider 
whitecaps to be any near-surface manifestation that results from breaking waves. Such 
manifestations include what could be called subjectively, “thick foam” (a lattice of bubbles 
floating on the sea surface), “thin foam” (a single layer of bubbles on the surface), and 
“subsurface bubbles” (air injected into the water during the breaking process). They do not 
include surface roughness normally associated with “sun glitter.” Historically, the whitecap 
coverage (or fractional coverage) was usually taken to be the fraction of the sea surface covered 
by thick foam.] Previous researchers have used a whitecap reflectance of 50% - 100% [Payne, 
1972; Maul and Gordon , 1975; Gordon and Jacobs , 1977] and have assumed it to be constant 
at all wavelengths. Laboratory measurements carried out by Whitlock et al. [1982] yielded a 
value of -55% in the visible part of the spectrum for thick foam in clear water, and showed that it 
diminishes with increasing wavelength beyond -0.8 gm, with 5% decrease at 0.85 nm, 10% at 
1.02 imi , and 50% at 1.65 nm. Koepke [1984] determined the effective reflectance of whitecaps 
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to be 22% in the visible, taking into consideration the decrease in reflectance and increase in area 
of an individual whitecap as it ages. He then combined this wavelength-independent estimate of 
the effective reflectance with the fractional area of the sea surface covered with whitecaps as a 
function of wind speed [Monahan and O’ Muircheartaigh, 1980] to obtain the overall reflectance 
of whitecaps for different wind speeds. Using Koepke's model for the reflectance as a function 
of wind speed, and the Whitlock et al. [1982] spectral variation of foam reflectance, Gordon and 
Wang, [1994b] carried out simulations that suggested that adequate atmospheric correction of 
ocean color sensors could be made at wind speeds -10-12 m/s even in the presence of large 
errors (-factor of 2) in estimating the contribution to the sensor radiance due to whitecaps. 

Recently, field measurements of the reflectance of foam generated by breaking waves in the 
surf zone [Frouin et al., 1996] indicated a larger spectral variation in reflectance than reported 
by Whitlock et al. [1982]. Frouin et al. [1996] found the reflectance decreased by 40% at 
0.87 nm, 50% at 1.02 (im, and 95% at 1.65 nm relative to the reflectance at 440 nm. The greater 
spectral variation in reflectance found in their field data, when compared to the laboratory foam 
measurements of Whitlock, is thought to be due to the strong absorption properties of water at 
the longer wavelengths. This absorption attenuates light reflected from submerged bubbles that 
are forced into the water column by the breaking waves. Using the spectral variation reported by 
Frouin et al. [1996], Gordon [1997] recomputed their influence on atmospheric correction and 
found whitecaps to be a far more serious problem than Gordon and Wang [1994b]. 

In addition to the spectral differences between foam measured in the laboratory and in the 
field, estimates of whitecap fractional coverage for different wind speeds are found to be quite 
noisy [. Blanchard , 1971; Monahan, 1971; Ross and Cardone, 1974; Wu , 1979; Toba and Chaen, 
1973; Monahan and O’ Muircheartaigh, 1980,1986]. Typically, fractional coverage (W) has been 
related to wind speed ( U) through an equation of the form W= a (A The values of a and p vary 
with geographic location, air and sea temperature as well as wind speed range. These 
relationships remain noisy, despite attempts to fine tune particular data sets, due to the dynamic 
nature and interdependence of the many parameters involved in whitecap formation. Whitecap 
coverage is primarily dependent on wind speed, but also on other factors such as fetch and 
duration [ Cardone , 1969], water temperature [Miyake and Abe, 1948], air temperature, stability 
of the lower atmosphere defined by the air/water temperature differential [Monahan and 
O’ Muircheartaigh, 1986], salinity [Monahan and Zietlow, 1969], and even surface tension 
variations due to the presence of organic films [Garrett, 1967]. A slightly different approach by 
Wu [1988] relates whitecap coverage to wind-friction velocity, being proportional to the square 
root of the wind stress which depends strongly on atmospheric stability conditions. Bortkovskii 
and Novak [1993] assessed whitecap and foam coverage with particular dependence on the sea 
surface temperature which effects sea water viscosity and wind friction. 

To date, the influence of whitecaps on the upwelling water-leaving radiance has been 
estimated by utilizing and combining the work of different researchers as described above. In 
this paper, a measurement system deployed from a ship is used to directly determine the 
augmented spectral reflectance resulting from the presence of breaking waves on the sea surface. 
Our measurement system provides spectral reflectance data from individual whitecaps — in 
various stages of growth and decay — and, by integrating over time, estimates of the augmented 
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spectral reflectance relevant to remote sensing can be determined. The principal reason we 
undertook the present investigation was to see if the Frouin et al. [1996] spectral measurements 
in the surf zone were applicable to real oceanic whitecaps. A secondary motivation was to 
validate the whitecap reflectance versus wind speed relationship proposed by Koepke [1984]. 

Instrument System 

The whitecap reflectance measurement system is described in detail elsewhere [Moore et al., 
1998]. The system consists of a 6-channel radiometer with narrow field of view and bands with 
a nominal 10 nm spectral width centered at 410, 440, 510, 550, 670 and 860 nm. It is held over 
the water surface by means of a boom extended from the bow of a ship providing an 
unobstructed view of the water surface. The field of view (FOV) of an individual channel on the 
sea surface is ~14 cm in diameter. Simultaneously, a deck unit with a cosine collector and 
matching spectral bands measures downwelling irradiance so that the reflectance of the water 
surface can be calculated. A TV camera is mounted next to the radiometer to provide a visual 
reference, and whenever the radiometric data are acquired the associated video frame is time and 
date marked and recorded onto video tape to assist in later analysis. Air and water temperature, 
wind speed and direction are measured simultaneously with the radiometric data at a rate of ~7 
times a second continuously for about 30 seconds (providing 200 contiguous samples), after 
which global positioning (GPS) data, universal time, and location are recorded. This acquisition 
sequence is repeated until a time determined by the operator. This radiometric rate allows 
several sample points of a large individual whitecap to be captured as well as providing an 
adequate data set over a time period of reasonably constant sea state and sky conditions. 

However, note that for a ship speed of 10 kts (5 m/s) the samples are -70 cm apart, so the 
radiometer will rarely view the same area of the ocean surface twice. 

Once the system has been installed on a ship, only periodic measurements of the dark current 
are necessary. This requires bringing the radiometer in from the boom and covering both the 
radiometer and the deck cell with light tight caps. This is necessary only when significant 
temperature changes are encountered such as when conditions change from early morning to 
noon, or from dark, overcast, rainy conditions to bright sunshine. The radiometer can also be 
angled up to 20 deg. from nadir in order to minimize the effects of specular sun and sky glitter. 

The purpose of the video camera was to remove events attributable to sun glitter from the 
analysis. Unfortunately, separating whitecaps from sun glitter was much more difficult than 
anticipated, and usually could not be done acccurately. Thus, most of the data reported here were 
acquired under overcast conditions (no sun glitter). 

Data Description 

Data used in this paper are from a -6,000 km transit through equatorial waters of the Pacific 
Ocean from Manzanillo, Mexico (19.03 deg. N, 104.20 deg. W) to Honolulu, Hawaii (21.20 deg. 
N, 157.55 deg. W). This 13-day cruise took place in November 1996 abroad the NOAA ship R/V 
Ka’imimoana. Conditions were uniform with air/water temperature and wind speed/direction 
changing only by a relatively small amount over time. Fetch was typically greater than 1000 km 
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and duration could be considered unlimited. Data were acquired far from land and for the most 
part, water temperature was -22-23C with air temperature -20-23C. Wind speeds ranged from 
~8 m/s to ~13 m/s and wind direction changed from an initial northerly direction at the beginning 
of the transect to an easterly and southerly direction as the Trade winds were crossed. Wind 
speed was constant within 1-2 m/s for days at a time. The sky conditions were mainly broken 
and scattered with some overcast periods. The data used in the analysis were collected under 
overcast conditions during 3 days of the transect. Reflectance data were acquired for wind 
speeds between -9-12 m/s. Stability of the lower atmosphere, defined by the air temperature / 
water temperature differential, varied from neutral to slightly unstable. 


Data Reduction Procedure 

The radiometric data collected were reduced to produce dark background-subtracted and 
radiometrically-corrected reflectance values, R, from the six up- and downwelling channels. This 
reflectance is defined as R(k ) = nL(k)/E(k), where L( A,) is the upward radiance leaving the sea 
surface, and E( X) is the downward irradiance falling onto the sea surface from the sun and sky. 

If the surface were lambertian (L independent of the viewing direction) then R would be the 
surface albedo. The radiometric correction process consists of three steps [ Moore et al., 1998] 
briefly reviewed here. 

First: ratios of upwelling radiance to downwelling irradiance are formed from the dark- 
subtracted signals for each channel. These six ratios are then multiplied by a radiometric 
absolute calibration factor relating the channels in the radiometer to those in the deck cell. This 
was established using a calibrated 1000W (FEL) quartz halogen lamp source and a calibrated 
Spectralon reflectance plaque. The deck cell was positioned 50 cm in front of the source and 
measurements taken. For the radiometer, the reflectance plaque replaced the deck cell, and the 
radiometer viewed the plaque at a 45 deg. angle. 

Second: a multiplicative factor was used to correct for spectral differences between the 
radiometer and deck cell and for differences between laboratory and field illumination 
conditions. 

Third: each channel of the deck cell was corrected for solid angle response for operation 
under natural illumination due to the slight deviation in angular response of the deck cell 
collector from a true cosine response. Correction factors for different angular distributions of the 
downwelling light field are required for the various sky conditions that may be encountered. A 
set of correction factors were established for sky conditions such as overcast, or clear with 
different solar zenith angles, to correct solid angle response. 


After applying these corrections we are left with a time series of spectral reflectance from 
the sea surface, examples of which are shown in Fig. 1. Figure la provides a 15 s record of the 
reflectance of two whitecaps passing within the radiometer FOV. It shows (1) that all spectral 
channels have essentially the same view of the whitecap, (2) that the reflectance of the whitecaps 
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throughout the visible is essentially independent of wavelength, and (3) that the whitecap 
reflectance at 860 nm is significantly lower than that in the visible, in agreement with the surf- 
zone observations of Frouin et al. [1996], Figure lb shows a longer record of the reflectance at 
410 nm. 

We wish to estimate the augmented reflectance [AR(A,)] associated with each whitecap event 
recorded in the time series. To automatically discriminate whitecaps and foam from the 
background reflectance we filtered each data set. The filter works in three steps on each 
wavelength individually. An example of the process steps is shown in Fig. 2. Because whitecaps 
are associated with rapid changes in the surface reflectance, the first step is to determine the data 
points for which the reflectance has changed by ±0.001 from the previous data point. For these 
data points, we believe there is a whitecap (or whitecap fragment) in the radiometer FOV. After 
this process is complete and candidate whitecaps have been determined, the background surface 
and subsurface reflectance must be eliminated. For each whitecap event we subtract the average 
of the first non- whitecap data point before and after the event. We are then left with individual 
cases of the enhanced reflectance (the augmented reflectance) due to the whitecap event. [This 
subtraction process inevitably leads to a few events with small negative augmented reflectance; 
however, their contribution to the final result will be negligible.] Finally we eliminate the events 
for which the augmented reflectance is less than 0.01. This last step, referred to as a threshold, is 
somewhat arbitrary; however, it results in fractional coverage estimates that are similar to those 
found by other investigators for similar wind speeds, e.g., Monahan [1971]. We will discuss the 
implications of this step below. Because the Fresnel surface reflectance is a much higher portion 
of the total reflectance for the red channels, the data are significantly noisier (especially at 860 
nm) and this procedure is much less reliable in the red than in the blue channels. This filtering 
procedure was validated by comparing the video record and the radiometric record. The filter 
(without the threshold test) passes candidate events that we subjectively consider to be 
whitecaps; however, it also can pass events that are only associated with changes in reflected 
skylight (due to the rough nature of the surface). These latter events are completely removed by 
the thresholding test. Thus, employing only steps 1 and 2, the contribution to the reflectance due 
to breaking waves will be overestimated. Using the threshold test (step 3) reduces this 
overestimation, and in fact, may underestimate their contribution somewhat. 

Within the context of ocean color remote sensing, the remote-sensing augmented reflectance 
(RSAR(A,)) is simply the increase in the water-leaving radiance (averaged over a single pixel, 
typically ~1 km 2 ) due to the presence of breaking waves on the surface. If the remote sensor had 
sufficient spatial resolution, one could in principle measure RSAR(A-) by subtracting the radiance 
measured in whitecap-free areas from the radiance averaged over a large number of pixels that 
include whitecaps and open water, after accounting for the transmission loss in the atmosphere. 
RSAR is then the average reflectance which the whitecaps have added to the surface reflectance. 
We attempt to effect the same determination at 410 nm by replacing the spatial average over -1 
km 2 by the time average of AR(410) over an ~14 cm diameter circle on the sea surface. The 
remote-sensing augmented reflectance at 410 nm was determined from the radiometer record by 
taking the sum of AR(410) for the reflectance “events” that passed these filters (Fig. 2) and 
dividing by the total number of records. Using the events that passed steps 1 and 2 of the filter 
provides an upper limit to RSAR(410) because examination of the video record shows that 
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events are passed that (subjectively) are not associated with whitecaps. These nonwhitecap 
events are removed (possibly along with some whitecap-associated events) by the threshold in 
step 3. Computation of RSAR(410) for the events that passed all three steps of the filter 
produces a lower limit to the augmented reflectance; however, as comparison with the video 
record indicated that the first-two steps of the filter often passed nonwhitecap events, we believe 
the lower limit is closer to the actual augmented reflectance. 

We can estimate the fractional coverage, W, by simply taking the ratio of the number of 
records containing events to the total number of records. This estimate will be higher than 
previous estimates (based on photography [Monahan, 1971; Koepke, 1984]) for two reasons: (1) 
steps 1 and 2 of the filter pass events that are not whitecaps; and (2) a small whitecap that is 
detected, but does not fill the radiometer FOV, will be included in the estimate of W as if it did 
fill the FOY. For the wind speeds observed here (9-12 m/s), W was -30% for events that passed 
steps 1 and 2, but fell to ~2-3% for events that passed the threshold as well. This latter range for 
W is more consistent with earlier estimates for these wind speeds. 

In Fig. 3, we present RSAR(410), as a function of wind speed, resulting from the analysis of 
1 1 data sets. Each of the 1 1 data points was for an average acquisition time of 45 minutes. 
Values of RSAR(410) that result from employing all three steps of the filter are indicated by the 

large triangles (A). To determine if the threshold in our filter was eliminating too much data, 
we also looked at RSAR(410) for which the threshold filter was set to 0. In other words we 
accepted all data for which the change in reflectance was greater than 0.001 from the previous 
point and for which the reflectance, after subtracting the background, was greater than zero. By 
comparison with the video records, this procedure appeared to include many low reflectance, 
nonwhitecap events, and should be a very conservative upper estimate of RSAR(410). Note that 
in this case the fractional coverage estimate increases to 30%, obviously too high, even 
considering the nature of the sampling. RSAR(410) for the zero threshold case is shown in Fig. 
3 as our upper limit (small filled triangles: ▼). 

To relate augmented reflectance to wind speed the power law relationship of Koepke [1984] 
has been modified from RSAR(410) = 6.49 x 10' 7 1/ 3 52 to RSAR(410) = 0.4 x 6.49 x 10' 7 I/ 3 ' 52 , 
where U is the wind speed in meters per second. This is the lower solid line in Fig. 3. Possible 
correlations between the augmented reflectance and other parameters such as water temperature, 
air temperature, stability of the lower atmosphere, humidity, salinity, barometric pressure, wave 
height, sea swell were investigated. None could be established mainly due to the small variation 
in the parameters during the transect. 

Because the filter scheme did not work reliably for the red wavelengths, another technique 
was used to determine the spectral variation of AR(A,). This was accomplished by passing the 
record for the wavelength X through the first two steps of the filter and comparing the result with 
the record for the 410 nm band passed through all three steps of the filter. If the filtered 410 nm 
band record indicated a whitecap, its augmented reflectance along with the augmented 
reflectance of the same event in the X band were recorded. Figure 4 compares the values of AR( 
A,) and AR(410) for individual events for one of the records that we examined. For each X 
(except 860 nm) the line on the figure is a 1:1 line. The noisy appearance in these figures is 
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partially due to the fact that each channel looks at a slightly different part of the surface. Clearly, 
for all channels (except 860 nm) the data fall near the 1:1 line. For 860 nm, the line shown is a 
least-squares fit (to an exponential) to the data. It clearly follows the 1:1 line reasonably well for 
low AR; however, for large AR, AR(860) is significantly less than AR(410). The depression of 
AR(860) at large values is similar to that noted by Frouin et al. [1996] for breaking surface and 
Moore et al. [1998] for ship wakes. 

Clearly, this analysis suggests that the reflectance of whitecaps at 860 nm decreases relative 
to that in the visible and that the fractional decrease is greater at larger AR(A-), i.e., AR(860)/AR( 
X) decreases as AR(A,) increases. The brightest whitecaps (here AR(A,)~0.3) show a reflectance 
decrease of -50% from the visible to 860 nm. Using the exponential relationship between 
AR(860) and AR(410) shown in Fig. 4, we can compute RSAR(860) for an entire record using 
AR(410) for each event to find AR(860), and averaging AR(860) over the record. The result is 
that for the range 9-12 m/s, RSAR(860)/RSAR(410)«0.8, with the ratio decreasing slightly as the 
wind speed increases. As the wind speed increases the occurrence of larger whitecaps increases. 
These larger whitecaps individually have the lowest ratio of AR(860nm)/AR(410nm), which 
lowers the overall ratio RSAR(860)/RSAR(410) for that wind speed. 

Concluding remarks 

We have determined the remote-sensing augmented reflectance in the open ocean for the 410 
- 670 nm range to be * 0.001 at a wind speed of 9 m/s, increasing to « 0.002 near 12 m/s (Fig. 3). 
These values are lower than previously derived values of whitecap optical influence (Fig. 3). 
Using the Koepke [1984] effective reflectance of 22% and the fractional coverage determined by 
wind speed for water temperature above 14C [Monahan and O’Muircheartaigh, 1980] yields 
values that are greater than 2.5 times our values from 9 to 12 m/s. Applying the fractional 
coverage relationship which takes into account atmospheric stability [Monahan and 
O’Muircheartaigh, 1986], our results are even lower than those for a stable atmosphere (Fig. 3). 
However, the Monahan and O’Muircheartaigh [1986] model does provide a significantly better 
fit to the data than the Koepke [1984] model. Also, the 860 nm value of RSAR, which will be 
dependent on the amount and depth of bubbles forced into the water column, is found to be 
approximately 20% lower than the RSAR in the visible for wind speeds 9-12 m/s in the open 
ocean. This is approximately half the depression observed by Frouin et al. [1996] in the surf 
zone, suggesting less violent wave breaking is manifest in the open ocean at these wind speeds. 

The goal of this study was to determine the augmented reflectance contribution from 
whitecaps as viewed from satellites and, in particular, its spectral dependence. Although most of 
our open ocean analysis utilized only data acquired during overcast conditions (to eliminate sun 
glint), some data for partly cloudy and clear skies were examined. This partly cloudy/clear sky 
data possessed an augmented reflectance spectrum similar to the overcast cases. Thus, we 
believe that our RSAR results should be reasonably representative of clear sky situations, and 
useful for estimating the whitecap contribution to the reflectance measured by ocean color 
sensors. 
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Notation 

AR(X) 
RSAR(X) 
Tair 
T water 

AT 

U 

w 


augmented reflectance at wavelength X. 
remote-sensing augmented reflectance at wavelength X. 
air temperature 
water surface temperature 
temperature differential (AT =T wa ter-Tair) 
wind speed in meters per second 
fractional coverage of whitecaps and foam 
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Figure 2. Sample record for Figure lb showing the operation 
of the filter at 410 nm. The solid line is the recorded surface 
reflectance. The dotted line is the background reflectance (re- 
flectance in the absence of whitecaps). 




RSAR(410) 



Figure 3. Remote-sensing augmented reflectance of whitecaps at 
410 nm. The large solid triangles are our measurements for white- 
cap events that have passed all three steps of the filter. The small 
solid triangles are our measurements for events that have passed 
steps 1 and 2 of the filter (no threshold). The open symbols 
are the Monahan [1971] fractional coverages multiplied by 0.22 
(the Koepke [1984] effective whitecap reflectance). The upper 
solid line is the Koepke [1984] whitecap reflectance model. The 
dashed lines are the Monahan and 0 ’Muircheartaigh [1986] model 
(also multiplied by 0.22) for a neutrally stable (AT = 0°C, lower 
dashed line) and an unstable (AT = 2°C, upper dashed line). 
The lower solid line is the Koepke [1984] whitecap reflectance 
model multiplied by 0.4. 








Figure 4. Comparison of AR( A) and Aft(410) for individual whitecap events. The solid lines on all 
graphs [except Ai?(860) vs. t1jR( 410)] are 1:1 lines. For the AR(860) vs. AR(410) graph (e) the 
curve is given by AR( 860) = 0.22 [l — exp[— 4.2 AR(410))]] . (a) AR(440) vs. AR(410); (b) AR(510) 
vs. AR(410); (c) AjR(550) vs. ,4i?(410); (d) Ai?(670) vs. Ai?(410); (e) AR{860) vs. AR(410). 








